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Abstract 



We analyze three classical field theories based on the wave equation: scalar field, electrodynamics 
and linearized gravity. We derive certain generating formula on a hyperboloid and on a null surface 
for them. The linearized Einstein equations are analyzed around null infinity. It is shown how the 
dynamics can be reduced to gauge invariant quanitities in a quasi-local way. The quasi-local gauge- 
invariant "density" of the hamiltonian is derived on the hyperboloid and on the future null infinity 
S + . The result gives a new interpretation of the Bondi mass loss formula Q, Q], We show 
also how to define angular momentum. 
^ ■ Starting from affine approach for Einstein equations we obtain variational formulae for Bondi-Sachs 

C*") ' tyP e metrics related with energy and angular momentum generators. The original van der Burg 

0] asymptotic hierarchy is revisited and the relations between linearized and asymptotic nonlinear 
situations are established. We discuss also supertranslations, Newman-Penrose charges and Janis 
' solutions. 

O ' 

1 Introduction 

o 

In the papers @, ||, [§) from the series "Gravitational waves in general relativity" Bondi, van der Burg, 
5-H . Metzner and Sachs have analyzed asymptotic behaviour of the gravitational field at null infinity. The 

energy in this regime so called Bondi mass was denned and the main property - loss of the energy was 
proved. See also discussion on p. 127 in the last paper JTo[ in this series. The energy at null infinity was 
also proposed by Trautman |Q and it will be called a Trautman-Bondi energy. 

We interprete their result from symplectic point of view and it is shown that the concept of Trautman- 
Bondi energy arises not only in gravity but can be also denned for other fields. In this case the TB energy 
can be treated formally as a "hamiltonian" and the loss of energy formula has a natural interpretation 
given by eq. ([l9]). We apply similar technique to define angular momentum. 

We introduce here the language of generating functions which simplifies enormously our calculations. 
This point of view on dynamics is due to W. M. Tulczyjew (see p^|). 

We start from an example of a scalar field for which we define TB energy as a "hamiltonian" on a 
hyperboloid. The motivation for concerning hyperboloids in gravitation one can find in ]l2| ], ]l5| | and ]lq |. 
In section 3 we give an example from electrodynamics. 

Next we prove analogous formulae for the linearized gravity. The result is formulated in a nice 
gauge-independent way. We show how the formula ( J19|) can be related with original Bondi-Sachs result 
- mass loss eqn (35) of || (cf. also eqn (4.16) in |3|, cqn (13) in Q and eqn (3.8) in jHj). Our 
result is an important gauge-independent generalization of this original mass loss equation. It shows the 
straightforward relation between the Weyl tensor on scri and the flux of the radiation energy through it. 
We show how to define angular momentum from this point of view. 

In section 8 we give covariant on a sphere formulation of the asymptotic equations from jij. We 
discuss several features of the theory like supertranslations, charges etc. and also the relations between 
linear and nonlinear theory. 
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2 Scalar field 



Consider a scalar field theory derived from a Lagrangian L = Li^tp^tp^), where ip^ :— d„tp>. The entire 
information about field dynamics may be encoded in equation 

6L(<p, <p„) = d^Stp) = (d^Scp + p^ . (I) 

The above generating formula is equivalent to the system of equations: 

n ,, dL „ dL 



Hamiltonian description of the theory is based on a chronological analysis, i. e. on a (3+l)-foliation 
of space-time. Treating separately time derivative and the space derivatives, 



sis, i. e. on a (3+1) 
, we rewrite (||) as 

5L= (ir6^y+d k (p k 5 l p) , (3) 
where we denoted ir := p . Integrating over a 3-dimensional space- volume V we obtain: 

5 L = / (irSp + irSip) + / p ± Sp = / (-hStp — (pSir + S(ttip)) + / p^Sip . (4) 
Jv Jv JdV Jv JdV 

where by p we denote normal part of the momentum p k . Hence, Legendre transformation between 7r 
and (p gives us: 

— SH(ip, 7r) = / (it8p — tpSn) + / p^Sip , (5) 

Jv JdV 

where 

Ti = f rnp — L . (6) 
Jv 

Equation (|J) is equivalent to Hamilton equations: 
6H . SH 

7T = -- r - ; <p>=— (7) 

op oir 

provided no boundary terms remain when integration by parts is performed. To kill these boundary terms 
we restrict ourselves to an infinitely dimensional functional space of initial data (tp, w) , which arc defined on 
V and fulfill the Dirichlct boundary conditions ip\av = f on its boundary. Imposing these conditions, we 
kill the boundary integral in (||), because Sp = within the space of fields fulfilling boundary conditions. 
This way formula (JsJ) becomes an infinitely dimensional Hamiltonian formula. Without any boundary 
conditions, the field dynamics in V can not be formulated in terms of any Hamiltonian system, because 
the evolution of initial data in V may be influenced by the field outside of V. 

Physically, a choice of boundary conditions corresponds to an insulation of a physical system composed 
of a portion of the field contained in V. The choice of Dirichlet conditions is not unique. Performing 
e. g. Legendre transformation between tp and p 1 - in the boundary term of (^]), we obtain: 

p^-Sip = 6 I p x (p - [ pSp 1 - . (8) 
av Jav Jav 

Hence, we have 

-S7l= [ {nS^p - pSn) - [ LpSp 1 - . (9) 
Jv Jav 

The new Hamiltonian 

H = H+ [ (10) 
Jav 



2 



generates formally the same partial differential equations governing the dynamics, but the evolution takes 
place in a different phase space. Indeed, to derive Hamiltonian equations (Q) from (^) we have now to 
kill Sp at the boundary. For this purpose we have to impose Neumann boundary condition p ± \ov = f ■ 
The space of fields fulfilling this condition becomes now our infinite dimensional phase space, different 
from the previous one. 

The difference between the above two dynamical systems is similar to the difference between the evo- 
lution of a thcrmodynamical system in two different regimes: in an adiabatic insulation and in a thermal 
bath (see (22)). As another example we may consider the dynamics of an elastic body: Dirichlet condi- 
tions mean controlling exactly the position of its surface, whereas Neumann conditions mean controlling 
only the forces applied to the surface. We see that the same field dynamics may lead to different Hamil- 
tonian systems according to the way we control the boundary behaviour of the field. Without imposing 
boundary conditions field dynamics can not be formulated in terms of a Hamiltonian system. 



2.1 Coordinates in Minkowski space 

We shall consider flat Minkowski metric of the following form in spherical coordinates: 

n^dy^dy" = -di 2 + dr 2 + r 2 (d6» 2 + sin 2 9d(f> 2 ) 

Minkowski space M has a natural structure of spherical foliation around null infinity, more precisely, the 
neighbourhood of S + looks like S 2 x M 2 . 

We shall use several coordinates on M 2 : s, t, r, p, uj, v, u. They are defined as follows: 

r = sinhcj = p^ 1 

t = s + coshw = s + /9 -1 yl + p 2 

u = t — r = s H , 

v = t + r = s + p _1 (VT+ p 2 + 1) . 



2.2 Scalar field on a hyperboloid 

We shall consider a scalar field <p in a flat Minkowski space M with the metric: 

n^dx^dx" = p- 2 (-p 2 ds 2 + J^M= + + d9 2 + sin 2 6dA (11) 

\ v/i + p 2 l + P J 

Let us fix a coordinate chart (x^) on M such that x 1 = 6, x 2 = <fi (spherical angles), x 3 — p and x° = s 

and let us denote by 7 ab a metric on a unit sphere (Jab dx A dx B := d9 2 + sin 2 #d0 2 ). 
We shall consider an initial value problem on a hyperboloid S: 

S s := {a; e M \ x° = s = const.} 

for our scalar field Lp with density of the lagrangian (corresponding to the wave equation) 

L := - i y/- det n^ jf v ip^ <p v = -^p^ 2 s'md 

We use the following convention for indices: greek indices p, v, . . . run from to 3; fc, Z, . . . are coor- 
dinates on a hyperboloid S s and run from 1 to 3; A, B, . . . are coordinates on S(s, p) and run from 1 to 
2, where S(s,p) := {x G S s | a; 3 = p = const.}. 



2, n2 (<Po) . 2(^3(^0 . ° AB 

1 + p 2 Vl + p 2 
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The formula (fL) can be written as follows: 



and 



5 L= (p», + / P 3 Sv 

JV JV JdV 



o 9L _ 2 / y? ^3 
P = — = P - 



3 9i -2 ■ n 1 

= = P sine 

OV3 V V 1 + P 



Let us observe that the integral / L is in general not convergent on E if we assume that ip = O(p) and 

Jv 

ip 3 — 0(1). We can "renormalize" L adding a full divergence: 



L := sin I 

2 



p\^f 



1 + p ; 



2 o^B 



V^+p 



L + \d (sing ll^ 2 ) - ^3 (sin0p-V 2 ) 



(12) 



where we have introduced a new field variable ip :— p (p which is natural close to null infinity. The 
generating formula takes the following form: 



S L= (tt°S^) + / tt 3 ^ 
Jv Jv ' JdV 



and 



dL 

dipo 

dL 

dip 3 



dL 



to 



V'3 



sin 8 



i + p 2 V^+p 2 
1 



= ^0 + P ^3 



AB 

— sin 7 



— - II 

It is easy to check that all terms are finite at null infinity provided ip — 0(1) and -03 = 0(1) 
From the above equations (in hamiltonian form) one can easily obtain wave equation: 

Dip = 

where the wave operator D is defined with respect to the metric 



rj ^dx^dx" := -p 2 ds 2 



2dsdp dp 2 



^TTp 2 i + p 2 

which is conformally related to the original flat metric t]^. 



d9 2 + sin 2 6»d0 2 = p 2 r] ^dx^dx" 



(13) 



(14) 
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Remark Let us observe that 



L = --^-det^Tf 1 '^ = L + \ (yp det^(ln p), u rj^ 2 



so we are not surprised that (13) holds. It can be easily checked that equation (|l^) is equivalent to the 
original wave equation: 



Ocp = 



(15) 



by the usual conformal transformation for the conformally invariant operator D + ii? because the scalar 
curvature R of the metric !) vanishes. 

If we want to have a closed hamiltonian system we have to assume that ip\ov — and then the energy 
will be conserved in time. But we would like to describe the situation with any data on S + ■ In this case 
we can define Bondi energy which is no longer conserved, formally we can treat it as a hamiltonian and 
formula (||) is useful as a definition of the Bondi energy together with its changes in time. In our case 
boundary condition / depends on time (see disscussion after formula (^)) and an interesting case for us is 
to compare the data with different boundary conditions. Although the energy defined on a hyperboloid 
is not a hamiltonian in a usual sense it plays an important role for the description of the radiation at null 
infinity. The metod is useful to the construction of the other generators of the Poincare group and will 
be applied for the angular momentum. 

Let us define "hamiltonian density" : 



H 



1 



and formally the following variational relation holds: 



S I H 
'v 



7T 3 5lp 



(16) 



dv 



here 7r :— ir°. 

Remark Relation between ( |l6| ) coming from L and (||) with respect to L can be described by following 
observations 

irSip — ipSn — p°5(p — Lp5p° 



it 3 Sip — p 3 Sip = -5 sinOpip 2 

So the formulae give the same hamiltonian because pip 2 vanishes on S + ■ 
We should express our hamiltonian as a functional of (Tr,ip): 



H := -sin( 
2 



TT^ 1 + P 2 



sinf 



\ D AB 

■03 +7 IpAlpB 



and the equations (0) are the following: 



= —^a+P 2 )+^vi+7 

sin 



n = (tiV 1 + p 2 )^ + [(1 +p 2 )sm<9V>3] i3 + (sin0 7 r/) B ),A 
and they correspond to the wave equation (fL3|) . 



(17) 
(18) 
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Although formally the formula ( [Tq ) looks similar to the usual hamiltonian system in our case there is 
no possibility to kill boundary term. Our "hamiltonian" is not conserved in time: 

-«9 (/#)=/ tt 3 V>= f sin6»(^) 2 (19) 

(9E=5(s,0) is odd oriented). 

Nevertheless this formal calculation is very useful. For example we can easily define angular momen- 
tum along the z-axis as a generator for the vector field -§t\ 



(n^5ip - tprfSn) = —8 / nip,,/, = —SJ Z 
Using equations of motion we can check that the angular momentum is not conserved in time: 



do J z = ~d ( / tt^ ) = / n 6 ^^ sinW,0 (20) 

'as Js(s,o) 



We will show in the sequel that the formulae (|19|), (|20|) can be written for linearized gravity and have 
also interpretation as the TB mass loss formula and angular momentum loss equation. 

Equations (|l^) and ( p0| ) are examples of the general formula which has the following form for a vector 



field X = 



do yj X U H + nX K tp k J = J nX^ip - n 6 X K ^ k - tt j X> (21) 

Remark The formula ( pl| ) defines usual hamiltonian system if the field X vanishes on the boundary. For 
example this situation take place for the so-called CMC (constant mean curvature) foliation which "ends" 
on the same sphere at S- 

We can take a boost generator along z-axis restricted to the hyperboloid Sq (with coordinate s = 0) 



X = -p^l + p 2 cos9d 3 - yjl + p 2 sm9d e 
and the formula ( ^l| ) takes the form: 

d ( f TrXVfe ) - / kX 3 iP - ir 3 X k ^ k = f sin 2 6^ e (22) 
Vis / J 9S Js(o,o) 



'S(O.O) 

Similarily the generator for linear momentum in z direction 
X = cos _ q q _ pi cos qq 3 — psinOde 

V 1 + p 2 

gives the loss formula: 



d P z = nX A il) - tt a X ^ k - 7r d X u V = - / sin 9 cos 9(ip) (23) 
ias Js(s,o) 



where P z := J X V H + TrX K ip k . 

The equations (p^|), (pp|), ( p^ ) and (^3|) express non- conservation law of the "hamiltonians" defined 
at null infinity. 

Let us formulate the following theorem: 

Theorem. If the TB mass is conserved than angular momentum is conserved too. 

This means that it is impossible to radiate away angular momentum without loss of mass. The proof is 
a simple consequance of (|l9|) and (pb]). If the TB mass is conserved than (from jlSj)) ip have to vanish 
on S and from (^oj) we get that angular momentum is conserved. 

We shall see in the sequel that this theorem also holds for Bondi-Sachs type metrics describing 
asymptotically flat solutions at null infinity for the full (nonlinear) Einstein equations. 



G 



2.3 Scalar field on a null cone 

We shall consider an initial value problem on a null surface N: 
N := jz € M \v = s + + \J\ + p 2 ) = constj 



(24) 



where we have introduced a null coordinate v :— s + p^ 1 (1 + \J\ + p 2 ) which plays a role of time in our 
analysis. Formally S + corresponds to the surface p — 0. Let us rewrite the Minkowski metric ( pd| ) using 
new coordinates v,u instead of s,p: 



r/^dx^dx" = p- 2 {-p 2 dv 2 - p 2 dvdu + d6 2 + sin 2 9d4> 2 ) 
The relation between coordinates (v,u) and (x°,x 3 ) used in the previous subsection is the following: 



v^x° + p- 1 (l + y / l + p 2 ), u=-2p-\ p = X 3 



d = d v , d 3 = 2p~ 2 d E - p- 2 1 



d v 



dx° = dv 



du , dx — —p du 



The density of the lagrangian takes the following form: 



L 



1 



detri^rj^ip^v = sin( 



1 oAB 



The formula ([[]) can be written as follows: 



S / L 
In 



N 



9JV 



and 



7T = — — = smdife 



= sin — = — sin 8 p 2 ips + 



Wo 



In this way we obtain density of the hamiltonian on a cone N 

„2 n AB 



H = tt v ^j v -£ = sm< 



(V*) 2 + j 7 VaV'b 



lim # = sin (9 (V^) 2 = sin6ty> 2 

p^0+ 



JV 



7rSip = / sin Oip^Sip 



A' 



sin OipSipdudOd^ 



(25) 



(26) 



We will show in the sequel that the above formulae exist in electrodynamics and linearized gravity. 
Similarity, the equation 



■Kip. 



A 



sin Oiprfip 



p^O 4 



AT 



sin Oipip^dudOdcf) 



(27) 



describes the flux of angular momentum through 
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2.4 ADM mass 



We have tried to treat separately hyperboloid and scri and we have learned that there is no possibility 
to get a nice hamiltonian system. Let N denotes in this subsection a "piece" of S + between E and i°. If 
we take the surface E U N together: 

-S[ I H+ I H) = (-kSip - ipSn) + / USip - ipSn) + / tt 3 # + / n^Sip (28) 

\JS JN J JS^ ' JN ^ ' J <9S JdN 

we will obtain hamiltonian system provided we can kill boundary term. This can be achieved assuming 
for example that 

lim ijj = 

u — y — oo 

which simply means that tj) is vanishing at spatial infinity. This usually hapcns for initial data on Cauchy 
surface t =const with compact support or vanishing sufficiently fast at spatial infinity. The following 
relations confirm our theorem: 

7r 3 | as = - sin0^> = 7r" 

<9£ = S(s, 0), dN = S(s, 0) U S(-oo, 0) 

— SrriADM = / I irSip — iPStt j + / sinOipSip (29) 

where m ADM ■= J^ uN H. 



2.4.1 one-parameter family of hamiltonian systems and their limit 

E Tj£ := {s = t, p> e} 



l + VTT 



e 2 e 



N r,e ■= < V = T + , == < P < 6 

e 1 + Vl + e 2 



Z T , e := it = t, < p < 



i + vm 



lim E r e = E T 
lim N T £ — N T C 

6^0 + 



lim I T c = i 



N T e := Iv = t-\ , h r < u < r 



i + VTT* 

t — r, r > 



E Ti£ U iV Tj£ U I T ,e is an explicit example of a one-parameter family of surfaces (with respect to r) and 
the hamiltonian related to this family is an ADM mass. On the other hand the hamiltonian system ( p8[ ) 
is a limit of these systems with respect to the second parameter e (e — > + ). This way we have certain 
"finite" procedure for the hamiltonian system d2l) at infinity. 
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2.5 Energy— momentum tensor 

Let us consider energy-momentum tensor for the scalar field tp: 
1 



-n 



where r\ := detrj^ and by 8^ v we have denoted Kronecker delta. In our case (for scalar field) symmetric 
and canonical energy momentum tensors are equal. 

For a Killing vector field we can integrate the equation 



and obtain as follows: 







do 



^{T\X V = - 



^T 3 V X V 



on 



It can be easily verified that for energy and angular momentum we have respectively 

/ v^T% = [ H 
Jn Jn 

/ V^V T °<j> = / TT*I>,<I> 
Jn Jn 

and the boundary terms can be expressed in terms of energy-momentum tensor 



T 3 p" 4 sin 



On 



T u vP ~ 2 sin 6d6d^ [ = 



on 



ip 2 sin 



A' 



T^p-^dpsinOdOdct) = I du[ -T v v p~ 2 sin 



N 



and ^P~ 2 T\ 



s+ 



Similarity for angular momentum: 



sin 



on 



T U & P~ 2 sin 



on 



9SC5 + 



4>4>,<l> sin 



f T v 4> p- A dpsm6d9d<t>= f • 
Jn Jn 



du -T v 



and ^P~ 2 T% 



i> , tp,< 



s+ 



Some of the equalities above hold for any energy-momentum tensor not only for the scalar field. Compare 
with (H^) and (f||-|27|) for scalar field. 

This calculation shows that quasi-local density of the energy on S + has two different interpretations. 
It is a boundary term which describes non-conservation of the "hamiltonian" on a hyperboloid £ or a 
density of a "hamiltonian" on <S + as a limit of N. More precisely, it is a density with respect to the 
parameter u but integrated over sphere. This is an example of an object which is local on Mi but 
non-local on S 2 . We call such objects quasi-local. It will be shown in the sequel that this concept of 
quasi-locality is useful in electrodynamics and gravitation. 
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3 Electrodynamics 



This section should convince the reader that the TB mass and angular momentum at null infinity can be 
described in classical electrodynamics in a similar way as the scalar field in previous section. 
Field equations for linear electrodynamics may be written as follows: 

SL = d ll (F , "8A v ) = d ll {T v ^)5A u + T v ^SA Vfl , (30) 

where A ViL := d^Ay and L is the Lagrangian density of the theory. The above formula (see J2j|) is a 
convenient way to write the Euler-Lagrange equations 

SL 



f) F v v — 

- SA„ 



(31) 



together with the relation between the electromagnetic field / M „ = A Vfl — A^ v and the electromagnetic 
induction density T v[1 describing the momenta canonically conjugate to the potential: 

(32) 



5A vli ' 

For the linear Maxwell theory the Lagrangian density is given by the standard formula: 



L = -K 

4 



and relation (^2) reduces in this case to T^ v := y/—r]r} tia r} l 'P f a p. 
Integrating (30) over V we obtain: 



(33) 



L = 



d {F k0 5A k ) + I T vi SA v = 



dv 



= f d a {T m 5A B + T 3t) SA 3 ) + f (T m 8A B + T° 3 SA ) . 

JV JdV 

We assume that the charge e defined by the surface integral 

S(s,p) 



(34) 



(35) 



vanishes. The situation with e ^ can be described similarily as in [B4| but we are interested in "wave 



degrees of freedom and we are going to show how the volume part of ( J34| ) can be reduced to the gauge- 
invariant quantities. 

Let a := p~ 2 A where A denotes the 2-dimensional Laplace-Beltrami operator on a sphere S(s,p) 
and one can easily check that the operator a does not depend on p and is equal to the Laplace-Beltrami 
operator on the unit sphere S(l). Operator a is invertible on the space of monopole-free functions 
(functions with a vanishing mean value on each S(s,p)). 

Let us denote by e AB the Levi-Civita antisymmetric tensor on a sphere S(s, p)). We can rewrite ( |34| ) 
in the following way provided that the electric charge e vanishes. 



[ L = f d [T" B , Ba - 1 Sp- 2 A B llB +^SA 3 +T 0B ,ce B c a - 1 6p- 2 (e AB A AllB )] + 
Jv Jv 



dV 



(^ B , Ba -Hp- 2 A B llB +^ 03 SA +^ B ,ce B o a - l S p - 2 (s AB A AllB )) 



(36) 



Here, by "||"we denote the 2-dimensional covariant derivative on each sphere S(s,p). Using identities 
d B !F m + d^J 730 = and d B T B3 + doJ- 03 = implied by the Maxwell equations and again integrating 
by parts we finally obtain: 



L 



do 



dV 



^ u a-^(aA 3 -(p-^|| B ), 3 ) 

||B,0 ) + 



{T 0B ^e BC )Br'5p-\e AB A A n B ) 



-1 s -2/AB 



^ 03 a-^(aA - p- 2 A B l{B , Q ) + (^^SBc^Sp-^e^A^) 



(37) 
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The quantities &A — p 2 A B \\ B0 and (aA 3 — (p 2 A B ^ B ), 3 ) are gauge invariant and it may be easily 
checked that 

sin0 (aA - p- 2 A B {{Bfi ) = p 2 T A Q \\A (= tt 3 ) 

and 

sin0 [slA 3 - (p- 2 A B llB ), 3 ] = p 2 T A 3 \\ A (= tt) 
Let us introduce the following gauge invariants: 
^ := jr 30 /sin6» 

*V := p-^^^sHA = - * F™l sm9 
tt := -p 2 T 3 A \\ A 

where 

Now we will show how the vacuum Maxwell equations 
= , M * = 

allow to introduce equations for gauge-invariants and the result is analogous to (jlj) and ([l8]) describing 
scalar field 



^=—{l + p 2 )+^Wl + P 2 (38) 
smt) 



7T = (7rVl + p2),3 + [(1 + P 2 ) skills] .3 + sin0a^ (39) 



^= ^(1 + P 2 ) + *0 3 V1 + P 2 (40) 
sm " 



*7r = (wyl + p 2 ),3 + [(1 + p 2 ) smO *ip 3 ] + sin 9a *tp (41) 
The proof of ( p8|) is based on the following observations: 

sin^.o = .F 3 % = -.F 3A ,a = -^(1 + P 2 ) + Wl + P 2 



sm ( 



T3A -rr0^4 

TT = p J- 3 n A = 



± + p 2 

and 

sin 0^3 =^ 3 °,3 
Similarily for (^) we have the following relations: 

- (*T 0A llB e A B ) fi (*F 3A \\BSA B ), 3 + (*T AB llBC e A c ). = 

*?° A = -p 2 e AB T 3B 
*J** = p 2 e AB T 0B 
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*J- — p £ 7-03 

77 _ -4 -r-03 
-^03 — — P -r 

which allow to get an equation 



,3 



which is equivalent to (|39|). 

For (*?/;, *7r) the proof is the same provided we apply the Hodge dual * for the variables and equations: 

(tt, — -> (*7r, — — > {it, if)) 

Now we will show how our variables appear in formula (|37|). Let us perform the Legendre transfor- 
mation in the volume V 

~ [A 3 - a^Or^V),,] = -<5 [^ 03 (A 3 - a-^p-^nB)^)] + 

+ [A3-a- 1 (p- 2 A B || S ), 3 ] 5T° 3 

and on the boundary dV 

T™5 (A a-V- a A fl || fl) o) - ^ [^ 03 (A„ - a-V 2 A fl || B , )] + 

-(A -a-V- 2 ,4 B || B ,o) «^ 03 . 



This way the formula (37) may be written as follows: 

8 [ [L- flb(^a _1 7r) - fliCV'a - V^ A o||a)] = ~ [ do [t^- 1 ^ + *7ra- 1 5 * V] + 



^oH^a- 1 ^-^ 111 ^^*^ ■ (42) 
Finally we have obtained the following variational principle: 

8 L = - I d (na~ 1 5ip + ma^S *ip) + / ir 3 a~ 1 8ip + *7r 3 a _1 5*V' , (43) 

Jv Jv J dV 

where the lagrangian L is defined by 

L = L-d (^n) - d 3 tyaT 1 ??\ u ) . (44) 
and boundary momenta: 

tt 3 := = p 2 F 0A llA - = sinO (^L= + P 2 ^) 

m 3 := -F 3A U B e AB = p 2 *Fo A \\a = sind \^L=+p 2 *^ J 



From lagrangian relation ([13|) we immediately obtain the hamiltonian one performing the Legendre 
transformation: 



5 H = — 7ra l 8tp — ipa 1 8it + wra 1 8 ~ 1 5*ir+ I n 
Jv Jv Jav 



3 a^ 1 8ifj + ^a" 1 ^ *t[) (45) 



where 

H := —TTa~ l ij) — «ra _1 *ip — L , (46) 
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is the density of the hamiltonian of the electromagnetic field on a hyperboloid. 

The value of J v H is equal to the amount of electromagnetic energy contained in a volume V and 
defined by the energy-momentum tensor. 

t\ = r x f Xv + \szr x f KX 

We are not surprised that the quantity H is related to T° 
[ Hd9d<t>= [ V^T° o d6»d0 

JS{s,p) JS{s,p) 

and to prove it we can use the following identity 

p~ 4 sin 9 wa -1 *ip — -kaT 1 ^ — <9 3 (ifjaT 1 p 2 F A \\a) = 

= T m Tn3, - T aA \ p~ 2 Fn A \ \a - T aA ^ E 'e ab^ 1 Fqa\\b£ AB 
The non-conservation law for the energy we can denote as follows: 

-d [H=[ sin (^a" 1 ^ + *0a _1 */>)=- / Ua' 1 ^ + mjm' 1 *^) sin6>d6»d0 (47) 
Js Jan ^ ' Js(s,a) v 1 

For angular momentum defined by 

we have a similar formula 

— da J z — do / 7ra _1 ?/;.0 + ma' 1 *^ > = — I sin# lipar 1 -^^ + npa^ 1 *^ a ) d#d</> (48) 
Js ' Js(s,o) v ' ' 

but the relation with symmetric energy-momentum tensor is not so obvious. 

Using the relations 

7T = -P _2 sin0/ 3 A || A V,3 = -p~ 2 £ AB hA\\B 

*tt = p _2 sin6»*/3 A || A *V,3 = -p~ 2 e AB *f 3A \\B 
we can express J z in terms of (ir,ip,*!r,*ip) as follows: 

J z ~ J if; ^a _1 7r : + sinOd^a^ 1 *^ + J *t[j ^a -1 *^ — sin Od^a^np = 

= — na^ 1 ^^ + *7ra~ 1 «/> ^ + / sin Oipd^a^ 1 *^ 

where 8a ■= sa b 8b and we have used the identity 

/ sinOip ^d^a^ 1 *tp = — / sin 6mj)d^a~ 1 'il) 3 
Js 2 Js 2 

The boundary term J as sin Oipd^a -1 *if) usually has to vanish if we want to interprete the integral 
is V~ r lT 4> as an angular momentum generator but in our case ip, *tp do not vanish on S + and we obtain 
in general two different definitions of angular momenta J z and J z . 
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Let us observe that 
A = (p- 2 A B || B )^-^*V 

and 



f (^ A %) A = f sin^a-V" 2 ( AS ||B) 6~ I sin^a" 1 *^ 
jv ' Jav Jay 



so the angular momentum J z is related rather to the canonical energy-momentum tensor with gauge 
A B ||s = than to the symmetric one. More precisely, the canonical energy-momentum density T M „ is 
related with symmetric one as follows: 

r% = T X »A X . V - ±5» U L = yf^jT^ + {T^A V ) X 

For angular momentum we obtain: 

/ r% = - / 7ra" + wa"V,^ - / sm Otp A a~ l p- 2 A B \\ B = J z + 

JS Js ' J OS 

+ / ijj^ar l p-' 2 A B \\ B sinedeA(j> 

Let us observe that if tp^ and #p $ are vanishing on S + then J z is well defined in terms of the canonical 
energy-momentum tensor density T ^ and conserved. 

3.1 Electrodynamics on a null surface 

We will show now how the formula ( p6| ) can be obtained in classical electrodynamics. Let us consider 
volume V C N where N has been already defined by (^4j). 

f j cvv 5A v dpd6d<j> = [ F vp a- x Sp- 2 A B \\ B d6d<p+ 

JV JdV 

Jv Jv 



and similarily as on £ 

aA p - (p- 2 A B llB ), p = -^^F vA u = tp = 2 P~*t* 



sinf 



where u = — - and d p = 2p t\. For dual degree of freedom the similar relation holds 



T vAllB e AB = *T vA \ A = 2sin^- 2 *V,u 
We obtain gauge-independent part + boundary term + full variation 



f F vu 5A v dpd6d(p= [ r ,p aT 1 8p- 2 A B \\ B d6d<p - 8 [ sm6ipa- 1 ip tU dud9d<p+ 

JV JdV Jv 

+ / sinfl (ip iU aT 1 6ip + *0 jU a -1 (S *ip) duAOdcp (49) 



This equality means that modulo boundary term and full variation we can reduce our form on S + and 
the final form is similar to ( p6| ) and posseses quasi-local character, (it = v + u so on the surface v = const. 
we can use coordinate u as well as u and this observation refers to the objects on N but not on M). 
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Now we will show how the flux of energy through S + is related with energy-momentum tensor, 
similarily to section 2.5. 

Tv 1 — 4/r-vn\2 , — 4/, rvn\2 , „„ rvA rvB 

v = 7}P (/ ) + 2^ W ) + 2 VAB -> f 



/ T%p~ 4 sin6»dpd6»d0 = - / sin0dud0d0 (V^a" 1 ^ + "^a -1 *V>,„) 

+~ / p 2 sin0dwd0d0(V> 2 + *0 2 ) 
4 Jv 



the last term vanishes on scri 



= / du [sin 0d0d0 (^,„a- V.« + ^a" 1 (50) 

The integral on a sphere in quadratic brackets represents quasi-local density of the flux of the energy 
through S + ■ The main difference comparing with scalar field is that here there is no possibility to work 
with local density because of the operator a~ x and only quasi-local object assigned to a sphere can be 
well defined. 



4 Linearized gravity on a hyperboloid 

We start from ADM formulation of the initial value problem for Einstein equations Q . In subsection 1 
we introduce the hyperboloidal slicing and in subsection 2 we consider an initial value problem for the 
linearized Einstein equations on it. In subsection 3 we discuss "charges" on the hyperboloid and in the 
next two sections we introduce invariants which describe reduced dynamics. In subsection 6 we derive 
"hamiltonian" in terms of gauge invariant quantities. 

4.1 Hyperboloidal conventions 

The flat Minkowski metric of the following form in spherical coordinates: 

•q^dy^dy" = -dt 2 + dr 2 + r 2 (d(9 2 + sin 2 0d(f> 2 ) (51) 

with r = smhuj, t = s + coshw already defined in section 2, can be expressed in coordinates s,uj well 
adopted to "hyperboloidal" slicing of Minkowski spacetime M: 

r] iiv df l dy v = -ds 2 - 2 sinh uj ds duj + duj 2 + sinh 2 uj(d9 2 + sin 2 Odcf 2 ) (52) 

In this section we use different coordinate u) instead of p used previously but at the end we return to 
p to compare the results for the scalar field and linearized gravity. Let us fix a coordinate chart (y M ) on 
M such that y 1 = 0, y 2 = <fi (spherical angles), y 3 = u> and y° = s. So we have 

E s := {y £M : y° = s}= (J S.(u) where S a (w) := {y G S s : y 3 = lu} (53) 

and S s is a three-dimensional hyperboloid, S s (ui) = S(s, sin ^ - ) and 9S S = £3(00) = S(s,0). 

We use the similar convention for indices (as for coordinates (2^)), namely: greek indices p, v, . . . run 
from to 3; k, I, . . . are coordinates on S and run from 1 to 3; A, B, . . . are coordinates on S(r) and run 
from 1 to 2. 
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Hyperboloid S has a very simple geometry. The induced Riemannian metric rjki on S in our coordi- 
nates takes the following form: 

r) k idy k dy l = duj 2 + sinh 2 Lu{d9 2 + sin 2 Odcf) (54) 

The hypersurface S is a constant curvature space and the three-dimensional curvature tensor of £ can 
be expressed by the metric: 

3 Rijki = VjkVa ~ VikVji (55) 



4.2 ADM formulation for linearized gravity on a hyperboloid 

Let (g k i,P kl ) be the Cauchy data for Einstein equations on a 3-dimensional hyperboloid E. This means 
that gki is a riemanian metric on S and P kl is a symmetric tensor density which we identify with the 
A.D.M. momentum ||, i.e. 



P H = Vdetg mn {g kl TrK - K kl ) 

where K k i is the second fundamental form (external curvature) of the imbedding of S into a spacetime 
M which is now curved. 

The 12 functions (g k i,P kl ) must fulfill 4 Gauss-Codazzi constraints: 



P/^ = 8n^det g mn T ifl n>* (56) 

(detg mn )H - P kl P kl + ^(P kl g k if = Wn{det g mn )T^n v (57) 

where T^ v is an energy momentum tensor of the matter, by 1Z we denote the (three-dimensional) scalar 
curvature of g k i, is a future timclikc four-vector normal to the hypersurface £ and the calculations 
have been made with respect to the three-metric g k \ ("|" denotes covariant derivative, indices are raised 
and lowered etc.). 

Einstein equations and the definition of the metric connection imply the first order (in time) differential 
equations for g k i and P kl (see j|] or @ P- 525) and contain the lapse function N and the shift vector N k 
as parameters: 

9kl = ^ (Pkl - \guPj + N k]l + N llk (58) 
where g := det g mn and P :— P kl g k i 

P kl = -Ny/g (n kl - ^g kl nj - ^ [p km p m l - ^PP kl ^j + (p kl N m ) lm + 

™ M ( r>kl — ' 



-9 



(p kl P kl - X -P 2 \ - N k [m P ml - N\ m P mk + Vff (n^ 1 - 5 H 7Vl" l | m ) + 



+ 87rN^T mn g km g ln (59) 

We want to consider an initial value problem for the linearized Einstein equations on the hyperboloidal 
slicing introduced in the previous section. For this purpose let us check first that on this slicing the ADM 
momentum P kl for background flat Minkowski spacetime on each hyperboloid S s is no longer trivial: 

P kl = -2^g kl (60) 

and 

9kldy k dy l = Vkldy k dy l = diu 2 + sinh 2 Lu{d9 2 + sin 2 9d<f> 2 ) (61) 
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where g kl is three-dimensional inverse of gki ■ 

Let us define the linearized variations (h k i,nj kl ) of the full nonlinear Cauchy data (gki,P kl ) around 
background data (|d0|), (|6l|): 



h ki := gki - Vki , vo Kl ■= P kl + 2Ar) M (62) 
where A := ^dct rjki (= sinh 2 a; sin #) 



We should now rewrite equations fl56|-p^) in a linearized form in terms of (h k i : vj kl ). Let us denote 
by w := rjki-uj kl and by h :— r\ kl hki- The vector constraint j56| ) can be linearized as follows: 

P l l \i~w l l \i-2Kh k \ k +Kh ]l (63) 

Let us stress that the symbol "|" has a different meaning on the left-hand side and on the right-hand 
side of the above formula. It denotes the covariant derivative with respect to the full nonlinear metric 
gki when applied to the P M but on the right hand side it means covariant derivative with respect to the 
background metric rjki ■ The scalar constraint ( |57| ) after linearization takes the form 

V9 n ~^= (P klp u - \{P kl 9ki) 2 ^j » A (h k \i - ftl fc ) |fe - 2vo (64) 

Linearized constraints for vacuum (!),„ = 0) have the following form: 

wi k \ k ~ 2Khi k \k + Ah\i = 0(= SnATi^ 1 ) (65) 

A (h kl \i - h} k \ - 2vj = 0(= le-KAT^n") (66) 

V / |fe 

The linearization of ([58]) leads to the equation: 

2N ( 1 \ 1 

hki = — l^fei - -jVki™ J + Kk\i + h i\k + 2Nr]ki{n + -h) - 2Nh H + 

—N m (h mk \i + h m i\k — h kl \ m ) (67) 

where N :— —J — = coshw, N 3 = 7703 = — sinhw, = j]oa = are lapse and shift for the background 

and n := is the linearized lapse. 



^kl 



= -NAh kl + Nw kl + N m vj k \ m + 2A (h klL + h llk ~ ?7 fe % m |m) + 



-A 



(Nn)\ kl ~ r] kl (Nn)\ m lm - (V lfc|i m + h ml \ k m - h kl \ m m - 



+ 



--N n 
2 



h kl\m + sh mi\k + sh mk\i _ rjki^lm + 2h mn ln )] (68) 



It is well known (see for example [p|) that linearized Einstein equations are invariant with respect to the 
"gauge" transformation: 

hfj, v + + £u;n (69) 

where £ M is a covector field, pseudoriemannian metric g^ v — rj^ + h^ v and ";" denotes four-dimensional 
covariant derivative with respect to the flat Minkowski metric ?7 M „. There is no (3+1) splitting of the 
gauge for hyperboloidal slicing similar to the situation described in [B. The (3+1) decomposition of the 
gauge acts on Cauchy data in the following way: 

A- 1 m kl A- X w kl + A^ 0|fe ' - N v kl ^ m m - 2Nfrj kl - N k ^ 1 - N l ^ k + 

+2r] kl N m £ Q \ m + 2^ L + 2^ k - 2r] kl £ m lm (70) 

h k i -> h kl +£ m + Z kll +2NrikiZ (71) 
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It can be easily checked that scalar constraint ( p6|) and vector constraint ( p5|) are invariant with respect to 
the gauge transformations (|70| ) and (f7l|). The Cauchy data (h k i, vj k[ ) and (h k i, W kl ) on E are equivalent 
to each other if they can be related by the gauge transformation ^ . The evolution of canonical variables 
vj k[ and hki given by equations ([37]), (^8|) is not unique unless the lapse function n and the shift vector 
h° k are specified. 

Let us define the "new momentum" p kl : 

p ki ._ w ki _ A ^ 2h ki _ ^klf^ ( p: =nj + Ah) 

This object can be also introduced in full nonlinear theory as P kl + 2^fgg kl and after linearization we 
obtain p kl : 

P kl + 2yfgg kl « p kl 

Let us first observe that the new momentum is trivial for flat Minkowski data. Secondly the symplectic 
structure is preserved: 

dP kl A dg kl - d (P kl + 2^gg kl ) A dg kl = -U 2 ^ = 
Moreover, the gauge transformation for p kl is simpler than for w kl : 

A -i p ki ^ A- 1 p kl +N^ kl -Nr 1 kl ^ m m -N k ^ l -N l ^ k + 2Ti kl N m ^ m (72) 
and the vector constraint has a familiar form: 

Wli = 0(= SttAT^) (73) 
We can also rewrite the dynamical equation ( |68| ) in terms of the new momentum: 

pki = N m p k ^ + N(r] kl p-3p kl )+A\(Nn)\ kl -r) k \Nn)\ m lm + 2Nnr) kl 

A 



+NA ( v kl h - 3h kl ) - —N (h mk \ l m + h ml \ k m - h kl \ m m - 
A. 



2^ 



il\k + h mk\l _ h kl\m + T] kl( h \m _ 2h mn \ n ) (74) 



We will show in the sequel that it is possible to define reduced dynamics in terms of invariants which 
is no longer sensitive on gauge conditions. The construction is analogous to the analysis given in 0. 

4.3 "Charges" on a hyperboloid 



The equation (|73|) allows to introduce "charges" related to the symmetries of the hyperboloid. There are 
six generators of the Lorentz group which are simultaneously Killing vectors on the hyperboloid S. Let 
us denote this Killing field by X h . It fulfills the following equation: 

X kll +X l]k = (75) 
Let V C E be a compact region in E. For example V := S s (r) and dV — S s (tq) U S s (ri). From 

re[r ,r!] 

(||) and @ we get: 

kl v / fnJvl v \ . I 



(Stt / AT^)=0= / p kl \ l X k = / (p fe 'X fc )| i= / p ik X k (76) 

JV JV JV JdV 

The equation ( [76] ) expresses the "Gauss" law for the charge "measured" by the flux integral. 

In particular for angular momentum when X — d/dip we can show the relation of this charge with 
dipole part of invariant y which we introduce in the sequel (subsection 4.4). 

16ns z := W<7rf y = -2 f p 3 = -2 / p 3 A {r 2 e AB cos(9)|| B = 

JdV JdV 
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= 2 1 r 2 p s A \ lB s AB cos 6 = / Aycosfl 
Jav Jav 

The time translation defines a mass charge as follows: 

(16tt f AT 0jU n") = = f N 
Jv Jv 



(77) 



K(h kl \i-h} k j fc -27r 



2JV fcP W |, 



2iV fc / z + A^A (V fc | fe - + A (A^/i fei - N l h) 

2N k p k3 + A (a/Vi 3 *^ - Nh 13 + N k h k3 - 7V 3 ^ 
and it can be related to monopole part of an invariant x (subsection 4.4). 



(78) 



16irp° 



av 



2N k p m + A (Nh 3k \ k - Nh 13 + N k h k3 - N 3 h"j 



A 



3V 



sinha; 



2 cosh 2 w/i 33 — coshw sinhwTJ, 3 — H 



2 sinli 
A" 



,33 



A 



av 



(79) 



Remark Traceless part of h k i and p kl has nice properties, gauge splits into 0-component (transversal 
to E) which acts on p kl and space components (tangent to E) which act on h k \ — ^r) k ih. But h and w 
remain nontrivial unless we impose gauge conditions. The most popular gauge condition, which allows 
to obtain scalar constraint as a full divergence (see below (|8C])), is to assume that -co = 0. Assuming this 
gauge we can define different "mass" charge as a surface integral coming from the scalar constraint ( |o6|) 
(but we obtain totally nonlocal object). More precisely, one can analyze the scalar constraint (|6(]) in the 
same way as (|76|) 



2 I w 
'v 



A 



Ah 



av 



(80) 



and there is no "Gauss" law for the "mass" defined by the surface integral of the right-hand side of ( p0[ ) 
unless we impose gauge condition zu = 0. This means that such definition of the mass charge measured 
by the flux integral at null infinity is not gauge invariant like ADM mass at spatial infinity. 



4.4 2+1 decomposition and reduction 

Now we introduce reduced gauge invariant data on E for the gravitational field similar to the invariants 
introduced in 0. For this purpose we use spherical foliation of E (see equations ([53]) and (|54j)). 

In this section we present mainly results without detailed proofs as in the section about electrody- 
namics. 

Let k := cothw. The gauge ( |7l| ) splits in the following way: 

h 33 -> &33 + 2&,3 + 2JV£° (81) 
h 3 A -» h SA + 6,a + U,3 ~ 2kCa (82) 
h AB -» h A B+U\\B + tB\\A + 2KriAB& + 2Nr)AB€ (83) 
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~N 3 (h 3AlB + h 3B \ A - h AB \ 3 ) = — pab ~ t^VabP + h 0A \\ B + (89) 



where by "||" we denote covariant derivative with respect to the two-metric r\ A B on S(r). Similarily the 
gauge (ff^) can be splitted as follows: 

A^V 3 -> A~V 3 + 7V<e° |33 - /V£ 0|m m = A- 1 P 33 -Nf^ A A -2Nnf' 3 (84) 

A -i p 3A ^ A-ip3 A + Ne l3A -N 3 e lA =A- i p 3A +Ne' 3A — m 

sinn oj 

A" V S -» A- V B + ^° |ylB - iV^ B C° |m m + 277 AB AT m ^l m = A- V B + (86) 

+m O\\AB _ NrjAB ^0,3 3 + C 0||C C + (K _ 27v 3 )e0j ^ 

It is also quite easy to rewrite (2+1) decomposition of (|S7|): 

fesa = ^ (paa - ^ + 2/i 03 |3 + 2/V(n + /i 33 ) - Af 3 /i 33 | 3 = (87) 

27V / 1 \ , 
= —(P^~2 P ) + 2h ° 3 ^ + 2NU + 2Nh33 ~ N 33 ' 3 
2N 

h 3A = —P3A + h 03lA + h QAl3 + 2Nh 3A -N 3 h 33lA = (88) 

2N o 
= ~P3A + h 03 , A + h 0AtS - 2nh 0A - N h 33>A 

2N ( 1 \ 

h AB = — I pab - -^VABPj +h A\B + h 0B \ A + 2Ni] AB n + 2Nh A B + 

2N ( 1 

x [ PAB - 2' 

+h B\\A + ^Kr]ABh 03 + 2Nrj AB (n + h 33 ) + 2Nh A B + 

-N 3 (h 3 A\\B + h 3B \\A - hAB,3) 

The vector constraint ( |73| ) can be splitted in the similar way: 

P3 k \k ^P3 3 ,3+P3 A \\A~ ^P AB VAB = (90) 
PA k \k =PA 3 ,3 +PA B \\B = P3A,3 + S A® \\B + 7} S \\A = (91) 

where S := p AB f] AB and S AB :— p AB — ^r] AB S. Similarily let us denote H := h AB r\ AB and xab ■— 
hAB ~~ ^VAbH. The invariants are defined as follows: 

x = 2cosh ujh + 2coshwsinhw/i nc + sinh uj\ \\ab — cosh u; sinh a; if, 3+ 

1/ . 2 sinh 2 u; ,0 . . 

-~(a+2)# p 33 (92) 

X = 2sinh 2 ujS AB \\ AB + 2 cosh uj sinh Lop 3A \\ A + ap 33 (93) 
y = 2A- 1 sinh 2 ujp 3A W B e AB (94) 
Y = A(a + 2)h 3A ^ B e AB - sinh 2 lo(A X C a\\cbS AB ), 3 (95) 
(2+1) decomposition of the scalar constraint ( |66| ) can be written in the form: 

A [h}\ - h kl \u?J +2m= [A(H, 3 - 2h 3A u - 2nh 33 + kH)] , 3 - 2K{h 33 + H) + 

+ 2(p 33 + S) - Hx AB \\ba + 2Kh 3 \ A ) + A (V 3 H^ + \h^ c ^J = (96) 
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The dynamical equations (f74|) take the following (2 + 1) form: 
A" V 3 = A- x iV3(p 33 ,3 - kS) - {Nn)^ A A - 2nNn, 3 + 

' [ h^ A A + H, 33 + 2kH, 3 - 2 K (2h 3A HA + h 33 , 3 ) - (2^^), 3 



±N 3 [H, 3 -2h 3A u ] 



A" V = A- x N 6 (p 6A , 3 + 2kjO + [(JVn), 3 - nNnf A + 



JV r 



3A ,AB. 



sinh cj 



|B,3 



2nh 



AB 



IIS 



IIS 



3A \\B 
1 \\B 



-N 3 h 33 ^ A 
2 



K^pab = A- 1 N 3 PABl3 +A- 1 N[r lAB { P 33 + S)+ PAB ] + 
+N(nu AB - T] AB n^ c c ) - Nr) AB (n, 33 + nn, 3 ) + 2rj AB N 3 n, 3 + 



+ - 



TV 3 
~~2 

N 
~~2 



hT A \\B + h* B \\A ~ VABh \\C + VAB(^H,a- h AA ]]A ~ h M ,a) - XA°,3??CS 



+ 



(X B,3VCA),3 + XAB U \\c~X a\\bc-X B\\AC + h 1 1 AB + -r)ABH\\ C + 



(97) 



(98) 



+VAB—^(h 33 + H) + t) AB (kH, 3 - 2 K h 3A u - 2K 2 h 33 ) 
Sinn u) 



—-2—XAB- 

smn u) 



r) AB (Kh 66 - -H, 3 ) + /i 3 A || fl + ft ) , 3 



We can check the reduced field equations for our invariants: 
x=-X + (iV 3 x),3 

X = N 3 X, 3 + AiVA s x - AJV 3 (x, 3 + 2kx) 

N N 3 
*=A Y + -^ 3 

Y = A(JV 3 A- X Y), 3 + AiVA s y - AiV 3 (y, 3 + 2 K y) 



(99) 

(100) 
(101) 
(102) 
(103) 



where A s is a Laplacian on a hyperboloid S. 

It can be easily verified that the invariants x and y fulfill usual d'Alambert equation (as a consequence 
of the above dynamical equations): 

□ x = 

□ y = o 

Let us notice that x and y are scalars on each sphere S s (r) with respect to the coordinates y A . 
For the scalar / on a sphere we can define "monopole" part mon(/) and a "dipole" part dip(/) as a 
corresponding component with respect to spherical harmonics on S 2 . Similarily "dipole" part of a vector 
v A corresponds to the dipole harmonics for the scalars and s v A u B . Let us denote by / mono-, 

dipole- free part of /. According to this decomposition we have: 



x = mon(x) + dip(x) + x 
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y = mon(y) + dip(y) + y 



Then mono-dipole part of the each scalar can be solved explicitly from the equations (100 103) and the 
solution has the form: 

4m 12k 



sinh u) sinh 2 ui 
12s 

y-y = . 2 

Sinn uj 
Let p := k than 

m = p = s = 

and 

k = p(s + cosh a;) + k 

Moreover am = 0, (a + 2)p = (a + 2)ko = (a + 2)s = which simply means that m is a monopole and 
ko, p, s are dipoles and they are constant on Mi. They correspond to the charges introduced in ||. Let 
us rewrite the solution in coordinates u, r which will be more useful in the sequel: 

, 4m+12p 12(ko+ptt) . . 

* = * + + ~2 ( 104 ) 

y = y + ^ (105) 

Let us remind the relation between spatial constant three-vectors in cartesian coordinates and dipole 
harmonics 

j lo zi p l zi s l z t 

k = , P = , s = 

r r r 

where zi are cartesian coordinates and j , p , s l are corresponding three- vectors representing our charges 
(see[|). 

4.5 Reduction of symplectic form on a hyperboloid 

We want to show the relation between the symplectic structure and the invariants introduced in the 
previous subsection. Let (p kl ,hu) and (s kl ,qki) denote two pairs of a Cauchy data on a hyperboloid. 

2 

The (2 + 1)- splitting of the tensor q^i gives the following components on a sphere: q:= r\ AB qAB-, Q33 ~ 

o 2 

scalars on S 2 , q^A ~ vector and <Iab '•= qAB — \f]AB 1 ~ symmetric traceless tensor on S 2 . Similarity we 

\ v r 

dipole part and the remainder in a natural way. 
The "monodipole" part we write separately: 



can decompose tensor density p kl . The quadratic form J v p kl qki can be decomposed into monopole part, 



mon( / p kl qki)= [ ^-9— p 33 mon(£) + / tan ^ ^ p 33 mon (g 33 )-|- 

Jv Jv 2 cosh uj J v A 



1 v 

If 2 

" tanhwp 33 mon(g) (106) 
dv 



2 



dip(/ P kl q kl )= I - \ 2 P 33 dip(0 - 2 / dip(sinh 2 V 4||S £A B )a- 1 (^|| B£ AB )- 
Jv Jv 2 cosh uj Jv 
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f / tan ^ ^ p^dip^ 33 ) + j- f tanho)p 33 dip(g)+ 
Jv A 2 Jay 

+ / dip ( — - — 5 h tanhwa _1 p 3A iu J ( -a 9 -2 sinh w cosha><7 3 ' 4 || J 4 

Jv V 2 cosh uj / \2 



v 2 cosh 2 w 

where invariant £ is defined as follows: 



(107) 



£ := 2 cosh 2 wq 33 + 2 cosh uj sinh wg 3 ' 4 || J 4 + sinh^ uj <l At> '\\ab — cosh uj sinh w 9, 3+ 



^3A 



2 , . XAB 



2 sinh w 



.33 



-2(a + 2)«-2" A - 
From vector constraints: 

sinhwp 33 , 3 + sinhwp 3A |u — cosh oj P= 



11^ 

,2,,; \u 



1 2 



(sinhW ,|a),3 + (sinh 2 uj P ab \\ab) + -aP= 



(108) 
(109) 
(110) 



(sinh 2 ^ 3A ll B eAB ), 3 + (sinh 2 ^ c Pa B \\bc) = 
we can partially reduce our form: 

J p kl Qkl= J p 33 q33 + 2p 3A q 3 A + ^PQ+P AB QAB = 

= / P 33 q33 - 2(sinhwp 3A || A )a" 1 (sinh^ 3 ' 4 || j 4) - 2(sinh ujp ZA ^ B e ab)^ 1 (sinh ujq 3A]]B e AB )+ 
Jv 

+ J i|?+2(sinh 2 we AC 'p A B || i 3 C )a- 1 (a + 2)- 1 (sinh 2 w£ AC 'g A B || BC )+ 

+2 / (sinh 2 wP AB || j4i j)a- 1 (a + 2)- 1 (sinh 2 w^ B || Ai 3) = 
Jv 

= / P 33 q33 - 2(sinhwp 3A || A )a- 1 (sinhtV' 4 || A ) - 2{sinhujp :iAllB eAB)a^ 1 {smhujq 3AllB e AB )+ 
Jv 

+ [ i(tanhcjp 33 ,3 + tanhwp 3A || A ) 1-2 [ (sinh 2 ujp 3A ^ B e AB ), 3a" 1 (a + 2)- 1 (sinh 2 uje ac q a B \\bc) + 
Jv 2 Jy 



(sinh ujp \\a),3 + -a(tanhwp , 3 + tanhwp \\ A ) 



^q -(a + 2)- 1 {smh 2 cjq AB llAB ) 



-i 

Jv 

= / tanhwp 33 
Jav 

-2 [ (smh 2 ujp 3A llA )a- 1 (a + 2)- 1 (smh 2 ujq AB llAB ) + 
JdV 

-2 I (smh 2 ujp 3A ^ B e AB ) a -\ a + 2)-\smh 2 u;e Ac q A B llB c) + 
Jav 



a- 1 (a + 2)- 1 (smh 2 ujq AB llAB ) = 



(sinh 2 ^ 3 ' 4 ll s eAB )a- 1 \q 3AnB s AB - (a + 2)- 1 (smh 2 ujs AC °q A B llBC ), 3 



+ 
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+ / p 

IV 



33 



g 3 3 + (a + 2)- 1 ( 



1 sinh 3 lo ° AB 



. 1 \\ab),3- -(tanhw«), 3 
coshu 11 2 



+ 



+ / tanhwp" 1 



- q +2sinhwcoshwa : (a + 2) ^sinh 2 a; 5 AS ||ab) ; 3 + 



-(a + 2) ^sinh 2 ^ Q AB \\ab) - 2a 1 (sinh lo cosh w<? 3A ||a) 

The volume term in the framebox is mono- dipole- free and corresponds to the invariants y, Y. The last 
two terms we can proceed further but let us first write a scalar constraint in two equivalent forms: 



.sinh lo o . R 

7— 1 AB \\ab ,3+ a + 2 

cosh LO ' 1 



1 2 

<733 - -(tanhr q), 3 



= (tanhw£), 3 + £+ 



2 

+ — (tanh 2 ws 33 — tanh oj sinh 2 lus 3A \\ A ) 



(111) 



ia(a + 2) q +2 sinh cj cosh w (sinh 2 lo <Z AB ||ab), 3 - a sinh 2 lo q AB \\ab+ 



— 2(a + 2)(coshw sinhwg ||a) = 



2 sinh lo 



= 2cosha>(sinhaj£), 3 — a£ — (as 33 + 2 sinh lo cosh ujs 3 \\ A ) 

In "radiation" part we get the following result: 



(112) 



.33 



. 1 .sinh lo ° 4 d , 1 , , 2 
933 + (a + 2) q AB UB ), 3 --(ta,nhu;q), 3 



tanhwp 3 ' 4 || j 4 



i q +2sinhwcoshwa *(a + 2) ^sinh 2 a; <1 AB \\ab), 3" 



-(a + 2) ^sinh 2 ^ 9' 4B || AB ) - 2a 1 (sinh w cosh w<7 3A ||a 

2 



(tanhu;£), 3 + £ + — (tanh 2 cjs 33 — tanhw sinh 2 c^s 3 " 4 ^) 



2 cosh w (sinh 3 — a£ 



2 sinh lo 



= f ^(a + 2)- 1 
Jv 

+ / tanhwp 3A || A a" 1 (a + 2)" 1 
Jv 

= / [tanhwp 33 + 2sinh 2 wa~V A ||A| (a + 2) _1 £+ 

JdV 

+ / [ap 33 - atanhwp 33 , 3 - 2sinhw(sinhwp 3 ' 4 || yl ), 3 - a tanh Lop 3A \\ A ] a _1 (a + 2) _1 £+ 
Jv 

+ f y [tanhcjp 33 (a + 2)~ 1 (tanh ws 33 - sinh 2 los 3A \\ A )+ 
Jv A- 

- sinh 2 Lop 3A | A a~ 1 (a + 2) " 1 (tanh was 33 + 2 sinh 2 los 3A | |a)] = 
= / [tanhwap 33 + 2sinh 2 wp 3A || A ]a" 1 (a + 2)"^+ 

JdV 



(as 33 + 2sinhwcoshws 3A || j4 ) 
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ap 33 + 2sinh 2 wP j4S || AB + 2sinhLJCOshL i ;p 3j4 || yl a 1 (a + 2)~ 1 £ 



+ J j [tanh 2 Lup 33 (a + 2)- 1 s 33 -2smh 4 Lup 3A llA a- 1 (si + 2)~ 1 s 3A llA \ + 

- [ ^tanhcjsmh 2 Lu(p 33 (a+2)- 1 s 3A llA +p 3A llA (a + 2)- 1 s 33 ) 
Jv A- 

and again framebox corresponds to the invariants but now x and X. 
Finally we get in volume integrand the gauge invariant part: 



p kl Qki = monodipole part in V + "radiation" part in V 



- / [tanhwap 33 + 2sinh 2 ^ 3A || j4 ]a _1 (a + 2) _1 ^+ 
Jav 

-2 f (sinh 2 ujp 3A ^ B e AB )a- x (a + 2)" 1 (sinh 2 ioe AC q 
Jav 



2 , ,^AC n B \ 
A \\BC) 



iav 

+ / tanhwp 1 
Jav 



:S.'S 



1 2 

2 



q-{a + 2)- 1 {s\x\h 2 uo°q AB \\ AB ) 



2 [ (sinh 2 wp 3A || A )a- 1 (a + 2)- 1 (sinh 2 w q AB \\ AB ) 
Jav 



where 



tanh lj 



p 33 mon(s 33 ) + 



onodipole part in V = / g — jJ 33 mon(£) + / 

Jv 2 cosh lu Jv 

f - \ 2 p 33 dip(g) - 2 / dip( S mh 2 Lop 3A W B e AB )a- 1 ( q3AllB e AB ) + 
Jv 2 cosh uj Jv 

[ tann ^ p^diprg 33 ) -f f dip ( — ^— = h tanhwa _1 p 3j4 H/i^ ( -a I -2 sinhw cosh wo 3A i U ^ 

7y A y ' Jv V2cosh 2 cj 11 y V 2 y 



"radiation" part in V 



ap 33 + 2 sinh z u P Ats \\ AB + 2 sinh u cosh wp' 



,2 , , „AB 



3A 



a- 1 (a + 2)- 1 C+ 



2 / (sinh 2 ^H s £AB )a- 1 \q 3AllB e AB - (a + 2)- 1 (sinh 2 uje ac '\ a B \\bc\\ 



4.6 Hyperboloid around S + 

Let us return to the coordinate p :— sin 1 h - . The metric on M takes the starting form ([ll]). It is convenient 
to introduce new canonical field variables similar to the variables for the scalar field and electrodynamics: 



* K :=p 1 x, tf„ := p l y 



n, 



x 



n. 



Y 



Equations of motion are the same for both degrees of freedom: 
1 



*T - *T,p = 



sin 6 



T = x,y 
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ri T - (n T V± + P 2 ),p = sinfl [a* T + ((1 + p 2 )^r,p ), P ] 



and they are similar to (|17|), (la) for the scalar field and fl38|), (39) for electrodynamics. 

The reduction of symplectic form from the previous section allows to formulate the hamiltonian 
relation in terms of new canonical variables: 



V / i> r a- 1 (a + 2)- 1 6n r - n T &- 1 (& + 2)- 1 6V r =6H- 
Jv 



T=x,y 



~t=x,y 



V / [n T y/l + P 2 + sin 0(1 + p 2 )*x,pl a -1 (a + 2)- 1 <51< T 
'av 1 J 



(113) 



where 



T=x,y 



tf T , P a- x (a + 2) 



sin 



- V] [ p 2 sin6»* T , p a x (a + 2) 1 *x )P - sin^ T (a + 2) 1 * T 
2 „ Jv 



T=x,y 



Similarily for angular momentum we propose the following expression: 
J z = J2 f n T a- 1 (a + 2)- 1 * T , 

T=x,y ^ s 

The non-conservation laws for the energy and angular momentum 
-d H= V / sin(9'i'Ta- 1 (a+2)- 1 'i<T 

•v- JS(0) 



(114) 



(115) 



T=x,y 



d J z = V / sin(9'i'Ta- 1 (a + 2)- 1 1'T, 
r=x, y J s(o) 



are similar to (|Ts]) , (^0|) and ( |47| ) . ( |48| ) . It should be also possible to formulate linear momentum P z in a 
similar way as (|23|) 



9 o p Z = J] / sin6»cos6'*Ta _1 (a + 2)" 1 'i'T 



T=x,y 

but this will be analyzed in a separate paper. It is obvious that all these formulae are quasi-local. 

4.7 Appendix — explicit formulae on a hyperboloid 

A = sinh 2 W sin0 N {m = -N m N {kl = -N k{l = N m N\ k k = -Nf k = ZN 

N = coshw N 3 = N 3 = - sinhw N A = ?7ab,3 = 2k?7ab 
77 AB , 3 = -2ki] AB £ab,3 = 2ke A b £ AB ,3 = -2ne AB n = cothw 

1 



A, 3 = 2kA kN, 3 = -kN 3 = N k, 3 = - 



sinh w 



r as — kijab r B3 = nS B r sc,3 = o 



2G 



RABCD — . : 2 {VACVBD - VADVBC) RAB = . : 2 VAB 

sinn u sinn uo 

£|33 = £,33 ^\3A = £,3A _ ^,A ^\AB = ^\\AB + ^AB ^3 
6|3 = 6,3 6|A = €a\\A - k£a £a|3 = ~ «&4 

€a\b = Ca\\b + nv AB ^3 

^33|3 = ^33,3 ^33|A = ^33,A — 2k/i 3 a 

h AB \ 3 — h AB , 3 + 2nh AB h AB \ 3 = 1iab,3 — ZkIiab h A B \ 3 = h A B,3 
h 3A l3 = h 3A , 3 + nh 3A h 3A \ 3 = h 3A , 3 - Kh 3A 

h 3 A\B = h 3 A\\B ~ kHab + K7]ABh 33 
h-AB\C = h A B\\C + KT]Ach 3 B + KrjBch 3 A 

5 Linearized gravity in null coordinates 

We are going to follow the idea from subsection 2.2 and apply it to linearized gravity. 
5.1 Minkowski metric in null coordinates 

Let us introduce null coordinates: u := t — r, v := r + t together with the index a corresponding to 
the coordinates (u,v). The spherical foliation is the same as previously and the coordinates on a sphere 
(x A ), (A = 1, 2), (x 1 = 8, x 2 = 4>) are the same. 

For convenience we need also some more denotations: p := r _1 = ^7^7, p, a = P 2e a where e u := |, 
e v := — i, rj ab e a e b = 1. We will also need e a := rj ab e b and we can check that e u = 1, e v = — 1, rj ab e a e b = 1. 

The explicit formulae for the components of Minkowski metric can be denoted as follows: 

—9 1 1 - 

VAB = P lAB, Vab = -^\E ab \, VaA = 

where E uu = = E vv and E uv = 1 = —E vu . Similarily the inverse metric has the following components 

AB 

^AB ^ p 2 1 ^ ^ab = _2\ E "b\ ^ ^A = Q 

where E uu = = E vv and E uv = 1 = —E vu . We shall also need derivatives 

■q AB , a = 2pe a V AB , VAB,a = -2pe a i]AB 
and finally the nonvanishing Christoffel symbols except T A bc are the following: 
T a AB = pe a VAB , T A aB = -pe a S A B 
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5.2 Riemann tensor in null coordinates 

We need to derive linearized Riemann tensor in null coordinates: 

2R a bcd — h a d.bc ~ hbd,ac + hbc ,ad ^ac.bd 

2R a bcD = haD,bc ~ hbD ,ac + hbc,aD ~ ^ac,6D + 

+peb (h a D,c + h c D,a — h ac ,D) - P£a {h b D,c + h cD ^ b — h hc , D ) 

ZRAbCd = h dA \\C,b + hbC\\A,d - Kd\\AC - h A c,bd+ 

+P e b {h d A\\C - h d C\\A - hAC.d) - P^d {hbC\\A - Ka\\C ~ hAC,b) + 
+pi]AC<^ a (hbd,a - h ad fi - h a b,d) - 2p 2 € b edhAC 

ZRABCd = h d A\\BC + h-BC\\A,d - h B d\\AC ~ ^AC\\B,d + 2pe d {hbC\\A ~ h A c\\ B ) + 

+pr)BC'£ a (haA.d - h d A,a + Kd.A + 2pe d h aA ) - pr)ACt a {Kba - h d B.a + Kd.B + 2pe d h aB ) 

2R a bCD = h a n\\C,b — hbD\\C,a + hbC\\D,a — h a C\\D,b + 

+2pe b (h aD \\ C - h aC \\ D ) + 2pe a (h bc \\D - h b D\\c) 

2RABCD = hAD\\BC + h B C\\AD - h B D\\AC - h-AC\\BD + 

+PVAC^ a {h B D,a - h aB \\D - h aD \\ B ) + prj BD € a (flAC,a - h aC \\ A - h aA \\c) + 
-pVBC£ a (hAD,a - h aA \\D ~ h aD \\ A ) - PVAD^ (^BC,a - h aB \\ c - h aC \\ B ) + 

+P 2 (h B DVAC + hAcVBD - hADVBC - h BC r]AD) + 2p 2 e a e b h ab (iiacVbd - VbcVad) 

5.3 Ricci tensor in null coordinates 

The Ricci tensor takes the following form: 

2R a b — h C b,ac + h a C \ c b — h a b' C c — h C c , a b + h a A,b^ A + hbA,a^ A — h a b^ A A — H ia b + 

+pt a Hj } + pe b H. a + 2pe c (h abtC - h ac ,b - hc,a) 

2R a B — h h B ^ ah - h aB c c + h a c iCB - h c CiaB + h a A \\ BA - h aB ^ A a + XB A \\A,a - 2^||s,a+ 
+pe a (2h b Bfi - h\ B ) - 2pe b h bB , a - 2p 2 e a e b h bB 

2RAB = (h a A\\B + h a B \\A) a ~ h a a \\AB - XAB' A a - 2pe a \AB,a + XA° \\CB + XB°' \\CA + 
1 



-XAB^ C C + VAB 



^H^c + H'\) + 2pe a {H ta - h a A ]]A ) + p 2 (2e a e b h ab - H) 
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5.4 Gauge in null coordinates 

The gauge transformation £ M 

splits in the following way: 

h a b * h a b + £ a ,b + (,b,a 

h a A > h a A + £,a,A + £,A,a + 2p€ a ^A 

hAB — > h AB + £a\\b + €b\\a - 2pVABe a £a 
and it would be also useful the following formulae 

XAB > XAB + £a\\B + £b\\A - VAB^ C \\C 

which are straightforward consequences of the previous one. 

5.5 Invariants 

Let us introduce the following gauge invariant quantities: 



y a := ( a + 2)h aAHB e AB -(p-z XA c llCB e AB ), a (116) 

y:= 2p- 2 {h bBU e AB ) ta E ab (117) 

x:= p- 2 X AB \\BA- + p~^H a - H + 2e a e b h ab -2p- 1 e a h a A u (118) 

x afc := a(a + 2)h ab - (a + 2) [{p- 2 h a A u ) , b + {p- 2 h b A \\ A ),a\ + [p- 2 (p~ 2 X AB \\ab),o\ m + 

+ [p- 2 ( P - 2 X AB UB ).b], a (H9) 



They fullfill the following equations: 

(p" 2 y a ),a = 2p- 2 y b , a E ab = (a + 2)y 

2E ab (p- 2 y), b + p- 2 y a = 

[p-\y a ,b - y b ,a)Y b + P~ 2 (a + 2)y a = 

(^ 1 y)' a a + P" 1 «y = o 

(p _1 x) l0 a = -p^ax 
p- 2 x ah , afc = a(a + 2)x 

71 ab *ab = (120) 
x ah = 2(p~ 2 x) :ab - r) ab (p~ 2 x)' c c 
if we assume vacuum equations = 0. 
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5.6 Reduction of symplectic form on S + 

Now we will show how linearized symplectic form on N, which was introduced by Kijowski in the full 
nonlinear theory, can be reduced to the invariants in "wave" part similarily to the hyperboloid case. 

We shall calculate this form in a convenient gauge but final result will be gauge-invariant. This way 
we shall prove that modulo boundary terms depending on the gauge the gauge-invariant part of the 
symplectic form can be obtained in the demanded shape in "wave" part. 

5.6.1 gauge conditions 

We would like to work with the following gauge conditions: 

Xab = 0, h a A \\ A = 
It is easy to verify that they are compatible for "wave" part: 

P- 2 X AB \\ab — P- 2 X AB \\ab + (a + 2)£V 

P-\a G \\cbs AB — P - 2 xa C \\cbs AB + { & +2)U\\bs AB 

K A \\A — > K A \\A + P~Va + (t A \\ A ),a 

More precisely, mono-dipole free part of £ Q and £4 is uniquely defined under these gauge conditions. 

5.6.2 partial reduction to extract gauge invariant part 

The linearized tt^ v has the following form 

n"" = -Ah' 11 ' + ^"(K + H)A 



and it can be simplified in our gauge. Let us observe that invariants (116) and (119) have simple 
form in our gauge in terms of h^ u . From (119) and ( |120| ) we obtain = —Ah uv — 0. Moreover 



r AB _ 

2' 

considerations we obtain: 



;r " = ^Arj AB h^ — Ax AB vanishes. Similarily ir \\a — because n Ab = —Ah Ab . Finally from the above 



[ 1L^5A^= [ K cd 5A? A - 2 ^ A W B e AB p- 2 S r l 5A a bA ^ 

Jsis.o) JS(s.o) 



'S(s,p) JS(s,p) 

One can show the following relation: 



/ K cd 5£ cd ~ I Ap 2 (p- 1 x afc ),„a- 2 (a + 2)- 2 <S(p- 1 x ah ) (121) 
Jv Jv 



where ~ denotes equality modulo boundary terms and full variation. 
Similarily one can prove 



f n bA ^e AB P ^ a -'S A\ AllB e AB ~ f Ap 2 y b a -\ a + 2)~ 2 6 [(p-^), b - {p'^] (122) 
Jv Jv 



5.6.3 full reduction to x,y 



We would like to obtain similar formula to (|26|). The "curl" part (122) reduces easily to the demanded 
form 

/ Ap 2 y b a - 1 ( a + 2)- 2 S[(p- 4 y v ), b -(p- 2 y b y v ] ~ f 2Ap 2 {p- 1 y) ^ ( a + 2)- 1 S(p- 1 y) 
Jv Jv 
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On the other hand the second part ( 121 ) can be rewritten in the following way 

/ Ap 2 (p- 1 x afc ), u a- 2 (a + 2)- 2 %- 1 x ah )~ f 2A / , 2 (p- 1 x), u a- 1 (a + 2)- 1 %- 1 x)- 
Jv Jv 



16Ap 2 



a- 2 (a + 2)- 2 ,5( i0 - 1 x) 



Let us observe that the last term vanishes on S + , more precisely (p 1 x) )tI = 0(p 2 ). Presented calculations 
should convince the reader that the following formula 

/ *rt£L~ I 2Ap 2 [(* ;c ) i „a- 1 (a+2)- 1 ( 5(v|/ ;c ) + (^) i „a- 1 (a + 2)- 1 ( 5(vI/ y )] 
Jn Jn 

holds and this is quasi-local form which is similar to (^6|) and (fl9|). 



6 Generating formula for Einstein equations 

There are different variational principles which may be used to derive Einstein equations. They may be 
classified as belonging to three basic approaches: 

1. The purely metric approach, where the variation is performed with respect to the metric tensor. As 
a Lagrangian one can use the second order Hilbert Lagrangian or the first order (gauge dependent) 
Lagrangian, quadratic with respect to the Christoffcl symbols. 

2. The metric-affine approach, based on the Palatini variational principle, where the variation is 
performed independently with respect to the metric tensor and to the connection. 

3. The purely affine approach where the variation is performed with respect to the connection. The 
metric tensor arises as a momentum canonically conjugate to the connection - see ]2^| . 

Each of these variational principles leads to the same Hamiltonian description of the theory. In this 
paper we use the Hilbert variational principle. At the end of this section we will show how the different 
variational principles converge to the same generating formula. Hence, the canonical structure, derived 
from this formula, does not depend upon the variational principle we begin with. 

The variation of the Hilbert Lagrangian 



may be calculated as follows: 



1 



5L = S[ — <T R»v ) = -j^Q^Sg^ + —y/\g\ g^SR^ 



where 



gfi u 



16tT 



16tt 



(123) 



(124) 



(125) 



We are going to prove that the last term in (124) is a boundary term (a complete divergence). For this 
purpose we denote: 



1 

167r" 



\9\ 9 



and 



A 



(126) 



(127) 
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(do not try to attribute any sophisticated geometric interpretation to A*; it is merely a frequently ap- 
pearing combination of the connection coefficients, which we introduce in order to simplify the derivation 
of the final formula). We have: 

®\ A nu = 9 x r^ u - d^r^ x = - r^ A r^„ + r^ CT r^ A (128) 

= Rfiu + A^ a A1 x — -A^A° a . 
Hence, we obtain an identity 

dx^SA^) = ^SR^ + ^dUl^-^A^AZ^j+idx^n^ (129) 
= ^ u 5R llv + {Vx^ u )5Al v . 
Due to metricity of V we have Va7t mi/ = 0. This way we obtain 

n^5R^ = d x {tt^SA^) = d K (*r K K») > ( 13 °) 
where we denote 



Inserting (130) into (124) we have: 



SL = -^-g^S gfil/ + dx(^SA^) . (132) 

We conclude that Euler-Lagrange equations Q^ v — are equivalent to the following generating formula, 
analogous to (|l|) in field theory: 

SL = d x (ir^SA^) (133) 

or, equivalently, 

SL = d K W K ^lu) ■ (134) 

This formula is a starting point for our derivation of canonical gravity. Let us observe, that it is valid not 
only in the present, purely metric, context but also in any variational formulation of General Relativity. 



For this purpose let us rewrite (132) without using a priori the metricity condition Va7t mi/ = 0. This way 
we obtain the following, universal formula: 

1 

167T' 



6L = -TT^Z^Sg^ - {V K *r K ) Kv + d - (^Ku) ( 135 ) 



It may be proved that in this form, the formula remains valid also in the metric-affine approach and in the 
purely-affine one. In metric-affine formulation, the vanishing of V x^^ is not automatic: it is a part of 



field equations. We see that, again, the entire field dynamics is equivalent to (134). Finally, in the purely 
affinc formulation of General Relativity the Einstein equations are satisfied "from the very beginning" 
whereas the metricity condition for the connection becomes the dynamical equation. We conclude that 



also in this case the entire information about the field dynamics is contained in generating formula (134). 

This formula, compared with (|l|) suggests that the role of field potentials in General Relativity should 
be rather played by the connection T, whereas the metric g should rather remain on the side of canonical 
momenta. This observation was the origin of the purely affine formulation of the theory. Also in the 
multisymplectic formulation (i. e. formulation in terms of Poincare-Cartan form - see |2^| ) the connection 
appears on the side of field configurations. We stress, however, that the results presented in this paper 
do not depend upon the choice of a variational formulation. 
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7 Metrics of Bondi— Sachs type 



In this section we shall consider the initial value problem for the curved space-time M with a metric of 
the form: 

g^dx^dx" = --^e 2/ W - 2e 20 dudr + r 2 j AB (dx A - U A du)(dx B - U B du) (136) 

on the null cone C = {x G M \ x° — u — const} (see j|) and boundary at null infinity dC — 

S(u, 0). We have the following non- vanishing components of the inverse metric g^ v : 

ff 33 = V e -2(3 

r 

g 03 = -e~ 2t3 

,,3A _ -2/3 TT A 



g^ = -c-^U 1 

9 AB = V 



where j AB is the inverse metric to ■jab- 

Let us define the "covector" Ub as follows: 

U B := 9baU A = r 2 lBA U A 

We have in our coordinate system the following non- vanishing components of the metric g^: 

goo = --c 2/3 + U A U A 
r 

903 = 

go a = -U A 

gAB = r 2 "/AB 



We also assume that 
\/detjAB = sin ( 



The metric (136) implies the following expressions for 167T7r'" y = y/^gg^ and A x fJ _ L , = T x ^ v — S^V 7 v ) a 
defined by (|J|) and (|l|7|). 



-g = e 2/3 r 2 sin ( 



167T7T 03 = -r 2 sin( 



AB 



16mr AB = e 2/3 sm9j 



167T7T 33 = rV sin ( 



167r7r 3yl = -r 2 sin 6U J 



167T7T 03 = -r 2 sin( 



167r7r AB =e^sin^ B 

^°33 = A° 3A = 
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^°03 = -0,3~- 

r 

A\ 3 = -- 
r 

A 3 3A = \e- 2P U B ,zg B A - \ (In sin ff) ^ 
A 3 o 3 = ~fi,a + (^) " uB fo -0- \^ P U A U A , S 

A 3 AB = \ e ~ 2P (y9AB ~ ^9AB,3 + U A \\ B + U B \\ A 

The following expression was proposed by Bondi to call the mass: 



m TB := ^- [ r-V (137) 

he 



8n 



Choose a (3+l)-foliation of space-time and integrate (133) over a 3-dimensional null- volume V <Z C = 
{x° = const}: 

S [ L= [ (n^SA%y+ f ^»5Al„. (138) 

JV JV JdV 

Similarily as in the case of electrodynamics, we use here adapted coordinates; this means that coordinate 
x 3 is constant on the boundary dV . Adapted coordinates simplify considerably derivation of the final 
formula. We stress, however, that all our results have an independent, geometric meaning. To rewrite 
them in a coordinate-independent form it is sufficient to replace "dots" by Lie derivatives Cx , where X 
is the vector field generating our one-parameter group of transformations which we are describing. In 
adapted coordinates X := Moreover, the upper index "3" has to be replaced everywhere by the 

sign _L, denoting the transversal component with respect to the world tube. This way our results have 
a coordinate-independent meaning as relations between well defined geometric objects and not just their 
specific components. 

Because the translation between these two notations is so simple, we have decided to use much 



simpler language, based on adapted coordinates. The volume part of the formula (138) can be simplified 
(or reduced) as follows: 



167T7r fcz ^; + 327r7r ofe (5A° fe + 167T7r 00 <Syl° 



= 327T7r 03 ( 5A o3 + IQ^ AB 8A° 



AB 



1 



--who (h ab \ s s {n AB ) + 



2r 4 sin ( 



(139) 



The last term in the above formula is a full variation of the quantity which logarithmically diverges 
when we try to integrate it, (3 — 0(r~ 2 ) and 2r 4 sin6> (jpi^j — 0(r _1 ). Removing of this term (- the 

renormalization of the symplectic form) corresponds to the renormalization of the lagrangian for scalar 
field (|12|). 
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On the other hand the boundary part in ( 138 ) can be rewritten as follows: 

lGnn^SA 3 ^ = 167T7r 33 (L4 3 3 3 + 327r7r 03 <L4 3 03 + 327r7r M <L4 3 3 A + l6irir AB 5A 3 AB = 

= 2sin0 (2V - r 2 U B llB ) 5/3 + sm6j AB 8U AllB -^sind (g AB - j9ab,3^\ 6~f AB + 



r 2 sin6e- 2f) U B t3 gBASU A - 6 



V 

-/ 

r 



rB /3,b 



(140) 



where we have denoted by " ||" a covariant derivative with respect to the two-metric g AB on dV. 
Inserting these results into (138) we obtain: 



KirrA / L = - j ^sinf 



(n AB \ 3 s( n AB ) 



JdV 



sm6e-^U B ^g BA 5U J 



[ 2sin0 (21/ - r 2 U B l[B ) 8(3 + )- sin0 (rV lAB , 3 - r 2 ^ AB - 2U A \\ B ) S-f AB + 

JdV 1 



l r 2 sin 9 


~2V + 








lav 







r 



(141) 



because 2-dimensional divergencies "d A f A " vanish when integrated over the boundary dV. 

It is convenient to introduce the following asymptotic variables (H AB ,tp AB ) related to asymptotic 
degrees of freedom: 



ip AB : = n AB 



Rab := --sin9 (rj AB ) 3 



From ( J139j ) we get the following relation: 



(IV 



~ sin 9 ( nAB ) 3 (n AB ) + 2r 4 sin 9 I A 



(142) 



On the other hand from we know that: 

16tt / n^C x Al v = f V^(V 3 X° - V°X 3 ) = 

Jv JdV 



r 2 sin I 



av 



y 

r 



- ) + 2-/3,3 - 2U B I3, B -2(3- e~ 2l3 U A U A 



(= 2r 2 sin 0^ 3 ) 



(143) 



where the last equality can be checked directly for the metric ( 136 ) and = Sq. From (142) and ( |l43| 
we obtain final formula: 



-S / 2V sin 9 + 
'av 







[ -sin9 


IV 


Jv 2 



av 



(n AB )6( nAB ) 3 -(n AB ) 3 6( n AB ) 

(rV^ AB>3 - r 2 j AB - 2U A]{B + r 2 e' 2f3 U c 3lcA U B ) 6-y AB + 



j 2r 2 sm9(^-U B llB + U A U A ^S(3-r 2 sm9e- 2l3 U A SU A (144) 
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Remark It seems to me that more natural "control mode" in the above formula corresponds to the control 
of the term {r 2 U A )^ than U A and it can be achieved by the following manipulation: 

-r 2 sm6e- 2l3 U A $Ui = - sm6 C - 2f3 U A S{r 2 U A )^ + S (r sin d&- 2fs U A XJ A ) + 



-2r sin 8e- 2f3 U A U A 5 (3 + - sin 9e- 2f3 U A U A S^ AB 
r 



If we pass to the limit the formula (144) takes the following form: 

- 16nSm TB = S [ 4Msin(9 = f tl AB 5ip AB - ip AB 5U AB -~ [ sm9r 2 ^ AB Sj AB (145) 
JdC Jc 2 J dc 

where V = r — 2M + 0(r~ l ) and the asymptotic conditions are given in and will be summarized in 
the next section. We can denote non-conservation law for the TB mass: 



1 ' • "„2- -,ab I 1 



167r9 m TB = / sm0r 2 j AB j AB [ = £ / sm9x AB , u X u ) (146) 

1 JdC \ Z JdC ' J 

where the last form in the brackets becomes clear when we learn about asymptotics presented in the next 
section. 

Similarily for angular momentum we get the answer from the superpotential proposed by Komar fl20f: 



16tt 



/ 7T^£ X A% = f V=9 (V 3 X° - V°X 3 ) 

JV JdV 



where now X — d/d(f>. 

The right-hand side can be expressed in terms of the Bondi-Sachs type metric: 

/ y/=g ( V 3 X° - V°X 3 ) = f r 4 sin 9e~ 2 ^ M U A — ► 16tt J z 

JdV JdV 

the limit is taken on S + and according to the asymptotics presented in the next section we obtain 

r i o bc 

16nJ z = - (67V + - X <f, B X ,, c )wa.66Sd4> 
Jdc 1 

But on the other hand 

167T / i^C x A% = [ n^A%^ = f T1 AB ^ AB A 

JV JV JV 

and 

16^0 / 7r"M° = f IW A % - H AB ^ AB = i / sm9r 2 j ABl AB ,^ 

JC JC ^ JdC 

We will show in the next section that non-conservation law for angular momentum agrees in terms of the 
asymptotics: 

f 1 o oAB 

16ttJ z = - -X AB:U X , sin0d0d0 (147) 
JdC 1 



3G 



7.1 Symplectic structure on scri 



Let us observe that we can use previous results (139) and (140) to reduce the form 

n^SA^ = n^5(A% + 2Al „) 

Let us also remind coordinate system which can be used to describe the situation in a similar way as in 
section 2.3 for scalar field and 3.1 for electrodynamics. (u,r) — > (v,u), u = — 2r v = u + 2r, d u = d v , 
d r = -2du + d v , du A dr = \du A dv and finally ir^6A v ^drd6d(f> = ^"SA^duded^. 
If we put 



Wtttv^SA ^ = ~ sm9( nAB \ 3 S ( n AB ) + 5 



2r 4 sin 6 



ft 



(148) 



and 



ltom^SA 3 ^ = 2sin0 (2V - r 2 U B \ {B ) 6ft + sm6j AB 5U AllB - ± sinO (g AB - j9ab,3 ) ^ AB 



! sm6e- 2l3 U B , 3 g BA 5U A -6 



2r 2 sin ( 



V 



V 

r 



+ -ft, 3 -ft-U B ft, B 



(149) 



assuming asymptotic behaviour on S + we obtain the following formula on the future null infinity: 



= - sinOr^^Sr^ 3 + 4<5(sin6»M) 



(150) 



Let N — [u.i,Uf] x S 2 C S + is a "finite piece" of S + . The relation with the TB mass is based on the 
following observations. First of all from (150) we obtain 

16tt f A v ^dudOd^ = -\ [ r 2 sin9^ AB ^ AB dudOdcj) + 2 f Msin 

Jn 2 2 J N J 9N 



and secondly 



16tt / n^C x A\ v 



N 



-g(V u X v - V"X") 



i)N 



2M sin 6d0d<t> 



ON 



where X = d , so finally 

-4 



Msin6d9d(f> = -- / r z sin 6*7^7 AB dud9d(j) 

dN 2 Jjv 



The left-hand side of the above formula represents the change of Bondi mass from initial state Ui to final 
state Uf (dN = {mj}xS' 2 U {ui} x S 2 ) but the right-hand side is a flux of the energy through N which 
is a piece of S + between initial and final state. 
Similarily for angular momentum we have 

16tt f ^"CxA\ v = 1&tx [ \^ v A\ v ^duded()> = -]- f ^sme^AB^duded^ 
Jn Jn 2 & Jn 



where now X = cL 
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8 Multipole structure of Bondi— van der Burg— Metzner— Sachs 
equations 



Let v = u + 2r than the metric (136) takes the following form: 

g^dx^dx" = (~ + r 2 lAB U A U B + e 2f3 ^j du 2 - e 2p dudv - 2r 2 lAB U B dudx A + 

+r 2 j AB dx A dx B (151) 

We shall rewrite formulae from van der Burg paper Q in a "spherically covariant" way. More precisely 
we denote: 

0. M i V are scalars 

1. pairs of functions U, W and N, P can be combined in two vectors U A and N A respectively: 

Ug = U 6 = U 

U$ = sin 2 OV* = W sin 6 

N e = N° = N 

iV = sin 2 0AT* = PsinO 

o 

2. pairs of functions c, d, C, H and D, K correspond to the symmetric traceless tensors Xab, Cab i Dab'- 

X g = -X = 2c 

X 4, = sin 2 6»Xe = 2dsin6> 

Similarily C e g — C, D e g = D etc. The reason of this notation arises in a natural way if we change 
the parameterization of the 2-dimensional metric jab- Let us remind that van der Burg in (p. 112) 
proposed the following parameterization: 

-f A Bdx A dx B = e 27 cosh(2(5)d6> 2 + 2 sinh(2<5) sin 6d9d<j) + e~ 27 cosh(2(5) sin 2 Ode/) 2 (152) 

which differs from original Sachs formulation by linear transformation of functions 7 and S (see H p. 
107). Next the used functions 7 and S are expanded as follows 

7 = c/r + (C - -c 3 - -cd 2 )r- 3 + Dr^ + 0(r- 5 ) 
6 2 

S = d/r + {H- -d 3 + -c 2 d)r- 3 + Kr~* + 0{r' 5 ) 
6 2 

Let us notice that there is no r~ 2 term which was analyzed in JTof and vanishing of this term is called 
"outgoing radiation condition" . 

We propose to change this parameterization in such a way that for original Bondi axially symmetric 
metric both formulations are the same and the main advantage of our change is that the expansion terms 
take a nice geometric form (mainly the term of order r~ 3 takes a nice form). 

o 

Let us fix the frame dO, sinOdcj) which is orthonormal with respect to the background metric Iab- 
The symmetric matrix (close to unity) 



e 27 cosh(2<5) sinh(2c5) 
sinh(2<5) e" 27 cosh(2(T) 

with determinant equal 1 can be also parameterized in a natural way by exponential mapping 

exp(aa x + ba z ) 



(153) 



38 



where a x and a z are Pauli matrices: 





' 1 o \ 




f 






v -1 J 


CT 2 = 




I 



The solution of the matrix equation 

e 27 cosh(2<5) sinh(2c5) \ , , . 

sinh(2«5) cosh(2£) ) = eXp(a ^ + 6 ^ 

leads to the relation between a, 6 and 7, <5 in the following form: 

arccosh(cosh(25) cosh(27)) 



(154) 



a = sinh(27) cosh(2(5) 
b = sinh(2<5) 



^/sinh 2 (2(5) + cosh 2 (25) sinh 2 (27) 
arccosh(cosh(2<5) cosh(27)) 



y / sinh 2 (2(5) + cosh 2 (25) sinh 2 (2 7 ) 
but the asymptotic relation for small 7, <5 is simpler: 
a = 2 7 + pS 2 + 0( l7 S) 5 

b = 2S+U 2 S + 0(j,S) 5 
o 

and it gives only correction in ?'~ 3 in our expansion. More precisely 

-a = c/r + (C- -(c 2 + d 2 )c)r- 3 + Dr~ A + 0(r- 5 ) 
2 6 

-b = d/r + (H~-(d 2 + c 2 )d)r- 3 + Kr~* + 0(r~ 5 ) 
2 6 

but now we can write the expansion in a matrix form: 

log7AB =Xab /r + (2C A b - ^ XcdX Xab) + 2D AB r~ 4 + 0(r~ 5 ) 



(155) 



and each term of the expansion is a traceless symmetric tensor on a sphere. The indices are raised with 

A R 

respect to the inverse 7 of the background metric (which is a standard unit sphere) . It is diagonal in 

00 o 

our coordinates lee= 1 and 7^0= sin 9. Metric connection of the Iab has the following non- vanishing 
components 

T 6 H = - sin 9 cos 9 , F% e = cot6l 

We are ready to show the asymptotic expansions for the rest of the quantities which appear in Bondi- 
Sachs type metric (136). They were introduced in Q (p. 114) but now we can rewrite them in a covariant 
way on S* 2 : 



U J 



1 oAB 

2r z ' 



2N J 



1 



loA BC 1 f CD 

2* bX ||c + — [XcdX 



U A := r 2 lAB U B = - 
i V _ 2M + N A M 1 



loB 2N A 1 

t;X a\\b H 1- -Tr- 

2 11 r lor 



,CD 



XcdX 



loAB o C 



1 oCDo 



£C \\bXa \\c + Yq X Xcd 
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1 1 o oAB 

f3 = ~32 ' r* XABX 
Basic equations, (eq. 13-15 in H): 

1 o AB \ AB 

M iU = --Xab,uX ,„ + 7 X \\AB,u (156) 
o 4 

4 114 1 oAB fo D EC\ 3 a A BC I CD A 

ZN A , u = -hf\ A - -e i^Xc \\de£ J --X b X \\c, u -^X ,„ X C \\d (157) 
(eq. 8-9 i 11-12 in j§): 

^ CD o 1 o CD o 

-4C*ab - - X XcdXab,« +- XabX X C d,u= N A]]b + N B{]A - 1 AB N c llc + 

o 1 o oC o p EG 

-MXab-^eacX bXe \\fgc (158) 
- AD AB = {a + 4)C AB - {Xa C ' N b )\\ c ~ (Xb C N a )\\c+ Iab (X ° N D \ C (159) 



Let us observe that mono-dipole but also quadrupole part of the right-hand side of the equation ( |159| 
vanishes. More precisely, from the relation 

rtAB\\C n AB C , nr ,AB 

<- CAB = t> ABC + \\AB 

we obtain that 

[( a + 4 ) C ' AB ]||AB = ( a + 6 ) (^11^) 



and similarily for Cab '■= £adC d B . Let us rewrite equation (15E) in the following way: 

-4Dab = {a + A)C AB - S A c B \\c (160) 

where 

Sabc -=Xac N b + Xbc N a - IabXcd N d 

It is easy to check that Sabc is traceless symmetric tensor (in each pair of indices) and the same holds 
for Sabc '■= £adS d B c- One can prove that S ABC ' \\cab and S ABC ' \\cab are orthogonal to the first three 
eigenvalue spherical harmonics (with I = 0,1,2). This way we get 10-dimensional space of quadrupole 
Newman-Penrose charges in Dab which are conserved Q). More precisely, quadrupole (and also 

'Vb and d u D AB lu 
asymptotics it may be not conserved (see PI). 
The linearized metric of Bondi-Sachs type 



mono-dipole) part of d u D \\ AB and d u D AB ' \\ AB have to vanish. However for the polyhomogeneous 



1 1 o oAB 

H=--^-X AB X =0 

2 r z 

Ku = l-^+2(3 + r 2 Iab U A U B - — + 
h U v = —/3 = 

t° a ° h 1 ° B Na 

h u A = -r 2 1abU b -rX AB U B = -- X A \\b + 

2 " r 

o o 1 2 ° ° 

Xab=tXab, h A B=rXAB+2 r ^ AB H = r Xab 
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and linearized asymptotics of invariants 



4M 6N A n A 



1 

"27 



, AB 



XabX 



lo c 
Y — --XA \\CB £ 



>AB 6 



,AB 



2^A\\B 



rx = AM , ^ x = rx.„ = X 



L4S,' 



3 o o^B l / CD 

^2 UcdX 



c 



ry = -Xa \\cb £ 



ry,u = 



c 



-,AB 



Xa 



\CB 



give an indication how to relate linearized theory with van der Burg asymptotics. This observation will 
be used in the sequel. 



8.1 Supertranslations 

Let us consider S + as a cartesian product S 2 x R 1 or rather trivial affine bundle over S 2 with typical fiber 
R 1 then the supertranslation corresponds to the null section of this affine bundle. On the other hand the 
boost transformation leads to the nontrivial scaling factor in a fiber and a conformal transformation on 
a base manifold S 2 (see || p. 111). 

Prolongation of the supertranslation from scri "to the center" in Bondi coordinates (metric ( |l36| )) 
leads to the following asymptotic relations (see also Q p. 119): 



x A + -a" A 



2r 2 



AB 

X a UR -2a llAB ai 



b + T A Bca""a 



\\BJ\C 



u = u — a 



Yr a ^ a W A + ^ 



, AB 



X ct\\AOt\\B — a' 



( Q 



B a< 



1 1 

r = r aa H 

2 2r 



oAB lo AB lo AB I 

X \\B a \\A + -jX a\\AB + -jX , u a\\A<X\\B - 2 a ®\\AB - Ot\\Att U " 



+ ^( aa ) 2 - ( aQ <*\\A 



Now we can check the transformation law for X and M: 

1 AB l AB l o AB 

M = M + -X\\ B ,u a \\A + 1 X „ "HAS + ^X, uu a nA a llB 

00 o 

Xab = Xab - 2a|| j4B + 7^saa 
Oa = Oa + a, aOo 
and finally we obtain that certain combination 



oAB o ab 

4M-X llAB =4M-X \\ab + a ( a + 2)a 

has a simple transformation law with respect to the supertranslations. Moreover, mono-dipole part of M 
is invariant with respect to the supertranslations. It corresponds to the mass and linear momentum at 
null infinity. This way we have proved the following: 
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Theorem 1. The energy-momentum 4- vector at null infinity is invariant with respect to the supertransla- 
tions. 

On the other hand angular momentum is not invariant with respect to the supertranslations but it 
transforms probably as follows: 

16tt7 z = 16ttJ z + / AMa 6 



Remark The supertranslation gauge freedom also exists in the linearized theory. The linearized part of 
the supertranslation corresponds to the gauge condition which preserves five components of the linearized 
metric: h ur , H, h rr , h r A- More precisely, it is a solution of the gauge conditions: 

£V = 0, r l|A A = (£V),r, (^Al\ B e AB ),r = 

II r 

We use here Minkowski background metric in the form: 

n^dx^dx" = -du 2 - 2dudr + r 2 7 ab dx A dx B 
and the solution of the gauge equations is the following: 

r " 2 

where a is any real mapping a : S 2 R and mono-dipole part of a corresponds to the usual translations 
in Minkowski space. The gauge transformation for tracclcss symmmctric tensor xab 

o 

XAB ► XAB ~ 2ra\\ AB + rl abSlu 

is similar to the nonlinear case. 

How to remove supertranslation gauge freedom? 

o ab 

Assume at time uq that AM— X \\ab — then stationary solution becomes simple stationary solution. 

o ab 

This procedure allows to treat X as invariant asymptotic degrees of freedom. 

Remark The Kerr-Newman metric in Bondi-Sachs coordinates can be asymptotically represented in such 

o 

a way that M = —Xab ■ 

8.2 Hierarchy of asymptotic solution on scri for scalar wave equation 

Let us rewrite wave equation in null coordinates (u, v) 

p- 1 (p- 1 i P y a a +a<p = (161) 



and suppose we are looking for a solution of the wave equation (161) as a series 

ip = (pip + ip 2 p 2 + f3P 3 + • • ■ (162) 

where each ip n is a func tion on scri, d v tp n — 0. 

If we put the series ( |l62[) into the wave equation ( |l6l| ) we obtain the following recursion: 

duVn+i = -7T-[a+ {n - i)n]<p n (163) 
2n 

Compare with equations 2, 3, 4 in 0. 

Remark The kern el of the operator [a + 1(1 + 1)] corresponds to the l-th spherical harmonics. The right- 
hand side of ( |l63| ) vanishes on the n— 1 spherical harmonics subspace. This means that the corresponding 
multipole in <p n +i does not depend on u. In particular for n = 3 we have quadrupole charge in the fourth 
order. The nonlinear counterpart of this object is called Newman- Penrose charge. We discuss some 
features related with NP charges in section 8.5. 



42 



8.3 Linear theory, asymptotic hierarchy, "charges" 

Let us first check that linearized theory can be obtained if we reject nonlinear terms in asymptotic 
hierarchy (157-159): 

4M =X \\ AB ,u 



3N A 



1 AB 



\\A 



-aM W A \\ B e =--&x AB \\ab 



-,AB 



-4C AB \\ab = (a + 2)N A \\A ~ 4C AB llAB , u = (a + 2)N A]]B e 
-4D AB \\ab = ( a + 6)C AB llAB -4D AB UB , U = (* + 6)C AB llAB 
x = 4M P + 6N A \\ AP 2 + 6C AB llABP 3 + 4D AB llAB p 4 + 0(p 5 ) 

AB 

y = X AB \\abP + 6N A{lB e p 2 + 6C AB llAB p 3 + 4D AB l{AB p A + 0(p 5 ) 
Full agreement with Jl63|) up to the 4-th order. 



oAB 

M = m + 3p + M, AM=X \\ab,u, rh = p = Q 



N A = -p» A u-k^ A -'e s llB + N" 
Let us notice that mono-dipole part of invariants: 



x = Amp + 12j lo xi ■ p :i + I2p l xi ■ p :i ■ [u 



y = 12s l x l -p i 



P 



x = 4(m + 3p) ■ p + 12(fc + p-u) ■ p 2 y = 12s ■ p 2 



is the same as (104) and (105). 



8.3.1 nonradiating solutions 



Suppose J s2 M = then from basic equation (156) we know 



= 



M = — 
s 2 8 



, AB 



Xab,u X 



and we get Xab,u= and finally also M = 0. Moreover equation (157) gives the following relation: 



,ab l 



CD 



3N A \\b s = -x" \\cd 

AB 

so the dipole part dip(iV j 4|| B e ) vanishes and this means that the angular momentum is conserved. 
This way we have proved the theorem formulated at the end of subsection 2.2, namely: 

Theorem 2. If the TB mass is conserved than angular momentum is conserved too. 

o 

The general solution of this type (namely M = =X AB . U ) will be called nonradiating solution and it has 
the following form: 



M = m + 3p + M 
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N A = -p^u - fcf - s ]]B + N A +~ (i I A B x CD HC DB Mil- 4 

4(7^ = -JV A ||b - N BU + 1 AB N C \\ C + (m + 3p + M.) Xab + \xabX CD \\cd 

^Dab = -{& + A)C AB + S A c B \\c 

Sabc =Xac N b + Xbc N a - IabXcd N d 

8.4 How to relate linearized theory with van der Burg equations, stationary 
solutions 

The "first order" asymptotics of the Bondi-Sachs type metric on S + is described by three functions 

o 

M, Xab- We shall try now to relate these data with boundary value on S of our invariants in 
the linearized theory in such a way that non-conservation laws for the mass and angular momentum are 

° 
similar in both cases. Suppose we know M, Xab at the moment u = uq and Xab,u on S (we need only 

in the neighbourhood of uo). We propose to perform supertranslation which is related with data at uq 

in such a way that 

Q AB 

(4M-Xiiab)(«o) = (164) 



Let us call the condition ( 164 ) supertranslation gauge at the moment uq. This way we have removed quasi- 

o ab 

locally supertranslation ambiguity at the moment uq. We stress that the relation AM— X\\ab = holds 

o oAB 

only at uq because QF(Xab,u) ■— 4 M u — X^ab u l& n °t vanishing in general however for the nonradiating 
solutions it may be fullfilled globally. This is the main difference between linearized theory where the 

condition AM— X\\ab = can be fullfilled globally and nonlinear data where we can only demand this 
condition to be fullfilled at one moment uo. Nevertheless these procedure allows us to relate nonlinear 
data at u$ with linearized theory namely 

AB o AB 

X^ab -> *y - X\\L,u - 

Now it is easy to verify analogy. The calculations from the previous sections devoted to the linearized 
gravity should convince the reader that in linear theory we can believe in the following equations: 

-16nd m TB = I sin0d0d0[* x a -1 (a + 2)- l ^ x + *j,a -1 (a + 2)" 1 * y ] (165) 

JS(s,0) 

-16nd Q J z = [ sin6»d6'd(/i[^a^ 1 (a + 2)- 1 * :Ei0 + * a a- 1 (a + 2)- 1 * 2/ .J (166) 

On the other hand the TB energy can be defined in terms of the asymptotics on S + 
16nm TB = / AM sin 0d9d(/> 



From (156) we get non-conservation law for the TB mass: 

f 1 f o oAB 

-167r5 m TB = - / AM sin0d0d</> = - / Xab,uX iU sin0d0d0 = 
Js 2 2 Js 2 
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sin edOdcj) 



l S 2 



o AB AB 

(X ll AB),ua- 1 (a + 2)- 1 (x llAB )^ + (XuB),na- 1 (a + 2)- 1 (xu B ),u 



(167) 



, AB 



Similarity between (|165| ) and ( |l67[) is obvious provided Xp B , n ^ X\\ab u ~ * 
Similarily angular momentum (around z-axis) can be defined as follows 



8ttJ 2 



S'2 



AB 

iVA||B £ cos sin 



where 



2Va := N A + X AB X ||c 



We have promised at the end of the previous section to show the relation (147) for angular momentum. 
The following sequence of equalities holds: 



16ttJ z = 



,BC 



S 2 

AB o 

X t u Xab, 



,BC 



1 oAB 



,BC 



6^ + -9„(X0sX 1 1 c) = / X^bX \\c,u + t;X ,« X<pa\\b ~^ X$b,uX \\ c 



where in the middle we have used equation (157). The last equality (in the above sequence) is a nontrivial 
identity and can be generally denoted: 



S 2 



-X xabx \\c + t;X Xab\\c-X - -X xabX \\c 



= -\ X AB {X c xab\\c + X c \\ b xca + X c \\ A xcb) 
1 Js 2 



(168) 



where now xab and xab are any symmetric traceless tensors on a unit sphere, X a 8a '■= and 

Xab, 4, = X c xab\\c + X c \\bXca + X c ' \\aXcb 
Another form of ( |168| ) can be transformed as follows: 

- / X BC (X A Xab\\c + X C XB A \\A - X a X bc\\a) = 
1 Js 2 

is equivalent to 

/ X BC X A {xbc\\a-Xba\\c) = I X BC X C XB J 
Js 2 Js 2 

and last equality holds for integrands 



\\A 



■BC vA /■ \ _ -BC vA - D _ Y A <y BC c c FD _ yA- FD 

X X \xbc\\a - Xba\\c) — X X eacXb \\d — X x ^ac^bfX \\d — X xafX \\d 



This way we have proved ( |168| ) and finally ( 147 ) which can be rewritten as follows 



— 1QttJ z = — 
2 



sin0d0d0 



oAB o 

X t u Xab,4>- 

, AB 



, AB 



(X|| AB ), u a- 1 (a + 2)- 1 (X|| AS ), + (x||i ; s ),ua- 1 (a + 2)- 1 (x||AB)^ 



The similarity with (|16q) is obvious provided that supertranslation ambiguity is removed. 
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Let us introduce the following objects in the full nonlinear asymptotics: 
4M = 4m+12p + Q + B x 

X AB \\ab = B x + Q X AB \\ab = By 

where m - monopole (mass), p - dipole (linear momentum), Q, Qo, B x , B y - mono-dipole free, 
represents supertranslation ambiguity: 

Qa = Qo - a(a + 2)a , Q = 



and equation (156) is equivalent to 

4m + I2p + Q = -\ X ABX AB = QF(B X , B y ) 

where QF is a quadratic quasi- local functional. Supertranslation gauge Q(u ) — Q allows to relate 
a t S + with Bx- The corresponding names m - mass, p - linear momentum, Q - "supertranslation 
charge" are obvious. This decomposition is chosen in a convenient manner for the situation of the so- 
called "sandwich- wave" . 

Suppose B x and B y have compact support on <S + (supp£?x C [ui,Uf] x S 2 C S + )- Let us also suppose that 
below m and upper uf our gravitating system is stationary. These two assumptions define "sandwich- 
wave" . 

When nonradiating solution becomes stationary? 

oAB 

From N A = we obtain p — M_ — X \\AB — 0, m, fco,s are not restricted but also X \\ab does 
not vanish. From Dab = we get Cab = (a + 4 : )~ 1 Sa c b\\c- Similarity Cab = gives m Xab — 

o - ~ AB oAB 

Na\\b + N B \\ A - Iab N^\\c or N A \\b £ =0 and X \\ab = (a + 2)N A \\ A - Let us call simple 

o 

stationary solution the situation when Xab— and M — m = const, described by van der Burg in static 
case (sec. 5 in Q). 

Remark The equations related with Newman-Penrose charge in static situation presented in Q at the 
end of page 119 can be denoted in our notation as follows: 

(a + 10)D AB = 15(MCab - N A N B + - 1 AB N c N c ) 
We have defined three categories of special solutions: 

simple stationary solutions C stationary solutions C nonradiating solutions 

and let us observe that the supertranslation gauge leads to the conclusion that every stationary solution 
in supertranslation gauge ( |l64j ) is simple. 

(nonradiating sol. in supertr. gauge) U (stationary sol.) = simple stationary sol. 

On the other hand in the case of the "sandwich wave" the supertranslation gauge at Ui and at Uf is not 
the same in general. The difference depends on 

Uf rllf 

Qdu= / QFdii 

so in general the initial and final states cannot be simple in the same Bondi coordinates. 
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8.5 Special solutions of asymptotic hierarchy, Newman— Penrose charges 



Equations ( |156| - |159| ) represent nonlinear analogue (up to the fourth order) of the hierarchy ( |163| ) for 
usual wave equation. 

We could define as a generalized NP charge any solution which starts in the n+l-th order from "multipole 
constant". More precisely, if tp S ker[a + (n — l)n] then from (163) <p„.+i.u = 0. Let us observe that if 
this charge vanishes we can derive "finite" Janis solution , which is obtained by "cutting the series" 
and derive hierarchy "upward" . More precisely 

If = Ifiip + lfi 2 p 2 + . . . + lfinP n 

ip n E kcr[a + n(n - 1)] => (p n+1 = , ip n = C(u)Y n _ 1 (6, (f>) 
2k — 2 

W - 1 = »(n-l)-(*-l)(*-2)** k ~ n 

and Yi is a spherical harmonics ([a+i(i+ l)]Yi — 0). In particular when (pi = then C(u) is a polynomial 
of degree n — 1. 

On the other hand if NP charge is not vanishing the solution ip can not be stationary. Moreover, monopole 
and dipole examples show that these solutions are singular (but on <S~). The monopole example is the 
following 

-i 2p 2 

ip = 2ip 2 pv =P2-z-> a( P2 = d u ip 2 = 

2 + pu 

similarily the dipole one 

4 3v — u 4 p 3 , 

^ = q^^7 \ = ~^ 3 q7^ ^2 a + 2(^3 = d n <p 3 = 

3 tr(v — ttj 3 (2 + pw) z 

and generaly 

ifn+i € kcr[a + n(n - 1)] => <p„+i = , ip n +i = CY n -i{0, 

but now (7 = and <p„ +2 = -^j<p n +i 7^ 0. 
For gravity we have: 

1. Dab = ==> Janis solution 

2. -Dyis - pure quadrupole NP charge solution. 

In Janis paper there are only linearized solutions. We shall try now to construct asymptotic quadrupole 
solution of nonlinear hierarchy (156 159|) . Let us assume that M_ and X CD \\cd are given quadrupoles 

(0 = (a + 6)M = (a + 6)x CI? ||cD)- 
M„ =Xab,«= 

o D CE 

4M = x(0,(t>) X C \\de£ =y{0,<t>) 
D AB = 

Sabc =Xac N b + Xbc N a - IabXcd N d 
Cab = (a + 4)~ 1 (S' y i C B||c) - (n A \\B +n B \\ A - ^ab n C \\ c 

N A = -p" A u - fcf - t B s llB + N A + V 
M = m + 3p + M 
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A 1 ° 

, := - e 




1 oAB 



CD 



\\cdb-M} ]A 



n 



N A \\B+N BllA - 1 AB N' 



He := (m + 3p + M) Xab + jXabX CD \\cd - 4(a + 4) H^sncJ 



This a special example of general nonradiating asymptotic solution denned by the condition that the TB 
mass is conserved. 



V - three-dimensional volume or function in Bondi-Sachs type metric 

L - lagrangian density 

S - hyperboloid 

ip - scalar field 

ip - rescaled scalar field 

5 - "variational" derivative 

<9 M - partial derivative 

T^v - symmetric energy-momentum tensor 

■q^v - flat Minkowski metric 

r] - det rj^y 

T M „ - canonical energy- momentum density 

S^u - Kronecker's delta 

p M - canonical field momenta 

7r - canonical momenta 

TL - hamiltonian, energy generator 

H - density of a hamiltonian or two-dim. trace of Kab 

,y" - coordinates on M 

t - time coordinate on M2 

r - radial coordinate on M2 

to - related radial coordinate on M2 , r = sinh u) 

p - "inverse" radial coordinate on M2, r — p^ 1 

s - "hyperboloidal time" coordinate on M2 , s = t — \/T + r 2 

u, v - null coordinates on M2, u = t — r, v = t + r 

u - coordinate on M2, u = — 2r 

0,4> - spherical coordinates on S 2 

d - exterior derivative 

/1, v, . . . - four- dimensional indices running 0, . . . , 3 

k,l, . . . - three-dimensional indices running 1, . . . , 3 

A, B, . . . - two-dimensional indices on a sphere 

a,b, . . . - two-dimensional "null" indices on M2 

D - d'Alambertian, wave operator 

EH - conformally related wave operator 

rj - conformally related metric 

7? - scalar curvature 



List of symbols 
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X - vector field 
i° - spatial infinity 
S - null infinity 
S + - future null infinity 
,S _ - past null infinity 
N - null surface "parallel" to 5 + or a piece of <S + 
rriADM - ADM mass 

S 2 - sphere parameterized by 6, <f> 
S(s, p) - sphere in M corresponding to coordinates s, p 
S s (lo) - sphere in M corresponding to coordinates s,u> 
S(s,0) - sphere on 5 + 
5(1) - unit sphere 

o 

Jab ~ metric on a unit sphere 

a - two-dimensional laplacian on a unit sphere 

oAB 

e ~ skew-symmetric tensor on a unit sphere, sin0e ^ = 1 
e AB - two-dimensional skew-symmetric tensor, r 2 sin de e ^ = 
| - two-dimensional covariant derivative on a sphere 
8a - dual of 8a, 8a = £a b 8b 
T^ v - electromagnetic induction density 
fnv - electromagnetic field 

- electromagnetic potential 
ip,*ip - gauge-invariant positions for electromagnetism 
7r,*7r - gauge-invariant momenta for electromagnetism 
J z - angular momentum 
gki - three-dimensional riemannian metric 
pki _ momentum 
K kl - extrinsic curvature 
R^vXcr - curvature tensor 
Rjiv - Ricci tensor 
r A M „ - Christoffel symbol 

n M - normal unit future directed vector 
1Z - three-dimensional scalar curvature 
h M - linearized metric 
vj kl - linearized momentum 
p kl - "new" linearized momentum 
g - det g kl 

A - volume element, A = r 2 sm9 
£ M - gauge in linearized gravity 
k - k — cothoj 
X, y, Y - invariants 

Xab ~ traceless part of Hab 
Sab - traceless part of pab 
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s 

H 

A E 

Jz 
Pz 

n«,n B 

T 

y« 

/3,\/,[/ A , 7 AB 

c 

mTB 

^AB 

Tl AB 
Cx 



P z ,P l ,P 

j l0 M 

M 
7,5 
U,W 
N,P 
N A 

c, C,D 

d, H, K 

o 

Xab , Cab , ^ab 
cr x , a z 
Sabc 
Sabc 
Cab 

XAB 

Dab 
mon(F) 
dip(F) 
£ 
F 
Y, 



trace of pab 
trace of Kab 

laplacian on a hyperboloid 

angular momentum generator 

linear momentum generator 

"asymptotic position" on a hyperboloid 

"asymptotic momenta" on a hyperboloid 

abstract index, T = x, y 

\E ab Ax a A dx b = du A dv 

invariant in null coordinates 

invariant in null coordinates 

parameters describing Bondi-Sachs type metric 

null cone or van der Burg asymptotics 

the Trautman-Bondi mass 

nonlinear asymptotic position on a null cone 

nonlinear asymptotic momenta on a null cone 

Lie derivative with respect to vector field X 

spin charge 

mass charge 

linear momentum charge 

static momentum charge (center of mass) 

Minkowski space or asymptotics of function V in van der Burg notation 
van der Burg parameterization of 7ab 
van der Burg parameterization of U A 
van der Burg parameterization of N A 
asymptotics of U A 

van der Burg notation for the asymptotics of 7 
van der Burg notation for the asymptotics of 8 

asymptotics of 7as 
Pauli matrices 



traceless symmetric tensor appearing in eq. ( 159 ) 

"dual" of Sabc, Sabc = £ a Sdbc 

"dual" of Cab, Cab =£a d C db 

"dual" of Xab, Xab =£a Xdb 

"dual" of Dab, Dab =£a C D C b 

monopole part of F 

dipole part of F 

mono-dipole free part of F 

supertranslation of F 

spherical harmonics with eigenvalue — 1(1 + 1) of the laplacian a 
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